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A  USERS  MANUAL 


FOR  COMPUTERIZATION 
OF  HYDROGEOLOGIC  DATA 
FOR  A  RIVER  ALLUVIUM 

By  James  W.  Naney,  Bill  B.  Barnes,  and  Douglas  C.  Kent1 

ABSTRACT 

A  system  for  rapid  storage  and  retrieval  of  data  from  a  multilayered 
aquifer  based  on  the  hydraulic  properties  of  the  aquifer  was  developed. 
The  manual  consists  of  the  data-retrieval  system,  control-card  setups  and 
specifications,  lists  of  variables,  formats  and  source  listings  for  data  stor- 
age and  retrieval,  and  graphical  illustrations  demonstrating  how  the  system 
works.  Data  for  the  examples  were  collected  in  the  alluvium  of  the 
Washita  River  near  Chickasha,  Okla.  These  data  were  processed  according 
to  the  manual  to  produce  a  lithologic  cross  section,  an  isopachous  map, 
a  specific  subdatum  map,  a  distributed  test-hole-data  map,  and  a  perme- 
ability-distribution map.  The  system  proved  capable  of  storing  and  re- 
trieving a  large  volume  of  layered  hydrogeologic  data,  and  it  will  have 
wide  geographic  application  where  hydrogeologic  studies  are  undertaken. 
The  system  may  also  be  adapted  to  other  types  of  materials  because  it 
uses  numerical  coding  exclusively  to  describe  the  hydrogeologic  data  to  be 
stored,  retrieved,  and  manipulated.  KEYWORDS:  computerization  of 
hydrogeologic  data,  ground  water,  ground-water  storage,  ground-water 
transmissibility,  hydrogeology,  hydrology,  lithology,  soil  permeability. 


INTRODUCTION 

A  layered  lithologic  system  created  the  need 
for  specialized  computer  programs  to  store, 
retrieve,  and  manipulate  hydrogeologic  data. 
Existing  programs  (7,  5,  7)-  would  store  data 
on  labeled  punched  cards  or  provide  for  the 
alphameric  coding  of  data  on  punched  cards 
and  store  and  retrieve  ground-water  hydro- 
graphs.  This  manual,  however,  has  been 
developed  for  the  selection  and  retrieval  of 
ground-water  data  based  on  the  hydraulic  prop- 
erties of  an  alluvial  aquifer.  It  documents  tech- 
niques described  in  an  earlier  paper  by  Kent  et 


'  Geologist  and  computer  specialist,  Agricultural  Re- 
search Service,  U.S.  Department  of  Agriculture, 
Chickasha,  Okla.  73018,  and  associate  professor  of 
geology,  Oklahoma  State  University,  Stillwater,  Okla. 
74074. 

-  Italic  numbers  in  parentheses  refer  to  items  in  "Lit- 
erature Cited"  preceding  appendix  A. 


al.  (3);  excerpts  from  that  paper  are  included 
here  for  convenience. 

The  retrieval  system  was  applied  to  define 
boundaries  and  internal  hydraulic  character- 
istics of  an  aquifer  in  the  Washita  River  al- 
luvium near  Chickasha,  Okla.  (fig.  1).  The 
alluvium  is  characterized  by  several  discon- 
tinuous layers  of  silty  clay,  sand,  and  gravel 
similar  to  those  found  in  stratified  alluvial  fans, 
as  well  as  some  stratified  basin  and  Coastal 
Plain  sediments.  Although  the  approach  used 
in  this  study  was  specifically  designed  for  al- 
luvial aquifer  systems,  as  shown  by  Naney  and 
Kent  (6),  it  can  be  used  for  other  similar 
aquifers. 

The  computer  codes  presented  in  this  manual 
are  written  in  FORTRAN  IV  for  the  IBM  360- 
65  computer.  Some  of  the  techniques  used  for 
processing  the  data  in  the  examples  were  de- 
veloped on  the  IBM  1130.  A  restriction  on  the 
word  length  for  coding  some  of  the  data  was 
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Figure  1. — Hydraulic-coefficient  ranges  in  the  Washita  River  alluvium. 


encountered  with  the  IBM  1130.  The  adaption 
of  these  programs  to  botli  computers,  however, 
has  expanded  the  availability  of  the  technique 
for  which  this  manual  is  written. 

Three  types  of  data  are  documented  and 
stored  in  this  data-handling  system:  (1)  gen- 
eral test-hole  information  (test-hole  data); 
(2)  stratification  characteristics  (layer  data) ; 
and  (3)  water-level  records  (water-level 
data).  The  test-hole  data  consist  of  the  well 
number,  location,  elevation  of  ground  surface, 
elevation  of  top  of  pipe  (well  casing),  the  num- 
ber of  layers,  and  the  method  of  drilling.  Layer 
data  include  layer  number  and  thickness,  meth- 
od of  analysis,  lithologic  type  and  color,  and 


approximate  range  of  the  hydraulic  coefficient 
for  each  layer.  The  water-level  data  include  well 
number,  record  type,  date,  time,  and  water-level 
elevation. 

DATA-RETRIEVAL  SYSTEM 

The  data-retrieval  system  (fig.  2)  is  based 
on  descriptions  of  the  lithology  in  terms  of  the 
type  and  color  of  material  and  of  the  hydraulic 
properties  associated  with  each  layer.  These 
descriptions  are  stored  by  test-hole  number, 
latitude  and  longitude,  and  layer  number  on 
a  disk  or  magnetic  tape.  General  test-hole  in- 
formation and  water-level  records  are  stored 
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Figure  3. — Flow  chart  for  retrieval  of  hydrogeologic  data. 


by  well  number  and  location  on  separate  disks 
or  magnetic  tapes.  Specific  information  related 
to  each  of  these  data  is  numerically  coded  and 
stored  by  variable  number.  Combinations  of 
these  stored  data  are  used  as  criteria  in  the 
selection  of  other  data  or  as  data  that  are  to 
be  shown  graphically. 

The  hydraulic  characteristics  of  a  multi- 
layered  aquifer  system  can  be  evaluated  by 
options  that  allow  selective  retrieval  of  data. 
Ten  options  for  data  selection  and  comparison 
are  specified  on  control  cards.  A  general  flow 
chart  of  the  retrieval  program  is  shown  in 
figure  3.  Figure  4  shows  how  each  specific 
option  is  selected.  Control-card  layouts  for 
retrieval  options  A  through  J  and  a  description 
of  card  types  A  through  N  appear  in  appendix 
A.  The  variable  name  and  columns  assigned 
to  each  code  are  shown  for  each  option.  A  list  of 
variable  names  with  descriptions  also  appears 
in  appendix  A.  In  addition  to  the  options  de- 
scribed, subroutine  packages  may  be  introduced 
to  be  called  by  specified  control  cards. 

Data-storage  format. — The  data-storage  pro- 
gram reads  data  from  punched  cards  in  the 
formats  shown  in  appendix  B  and  stores  them 
on  magnetic  tape.  Twenty-one  variables  are 
stored  for  each  test  hole,  including  variables 
that  describe  the  entire  well,  such  as  the  type 
of  well  casing  or  the  elevation  of  the  top  of 
the  pipe.  These  variables  are  presented  in  table 
1. 

Twenty  lithologic  zones  or  layers  are  allowed 
for  each  test  hole.  Currently,  no  more  than  14 
distinguishable  layers  have  been  stored  for 
any  one  test  hole.  There  are  eight  variables 
describing  the  lithologic  characteristics  of  each 
layer,  such  as  layer  thickness  or  the  type  of 
material.  These  variables  (table  1)  are  read 
from  punched  cards  according  to  the  format 
shown  in  appendix  B  and  stored  on  magnetic 
tape. 

Ground-water  levels  are  stored  on  magnetic 
tape  for  use  in  conjunction  with  the  test-hole 
and  lithologic-layer  data.  The  format  used  on 
punched  cards  to  store  the  data  is  shown  in 
appendix  B.  Storage  and  retrieval  programs 
for  test-hole,  layer,  and  water-level  data,  along 
with  comments  and  execution  cards,  appear  in 
appendix  C.  The  three  basic  graphical  forms 
produced  are  test-hole  cross  sections,  lithology- 
distribution  mans,  and  isopachous  maps.  A 
separate  retrieval  program  (not  shown  in  fig- 


ure 2)  provides  data  for  maps  that  show  the 
distribution  and  change  in  level  of  the  ground- 
water surface. 

Coding  and  storage  of  lithologic  data. — The 
lithologic  type  and  color  are  described  with 
unique  combinations  of  numbers  (table  2).  Only 
five  numbers  representing  grain  size  are  used  in 
the  code  to  describe  the  lithologic  type  because 

Table  1. — Data-processing  code  for  tape  stor- 
age, entire  well 


Variable  Term  Code  form1 


1  Well  number  or  test-hole 

number    XXXX 

2  Watershed  number    XXX 

3  Location   (latitude:  degrees, 

minutes,  seconds)    XX  XX  XX 

4  Location   (longitude:  degrees, 

minutes,  seconds)    XX  XX  XX 

5  Top  of  pipe  elevation    XXXX.XX 

6  Ground-surface  elevation    XXXX.XX 

7  Bedrock  elevation    XXXX.XX 

8  Total  depth  of  test-hole  log  (ft)  .  .  XXX 

9  Distance   from   logged  well  to 

observation  well  (ft)    XXXXXX 

10  Total  depth  of  observation  well 

■     (ft)    XXX 

11  Drilling  method    X 

Core    1 

Rotary    2 

Split  spoon    3 

Auger   4 

12  Casing  type    X 

Metal    1 

Plastic    2 

13  Casing  length  (ft)    XXX 

14  Casing  diameter  (in)    XX. XX 

15  Screen    XXX. XX 

Screen  length  (ft)    XX.XX 

Screen  diameter  (in)    XXX. XX 

16  Pump  type    X 

Submersible    1 

Jet    2 

Turbine    3 

Centrifugal    4 

17  Pump  use    X 

Domestic    1 

Irrigation    2 

Experimental    3 

18  Pump  capacity  (gal/min)    XXXX 

19  Number  of  stages    XX 

20  Pump  diameter  (in)    XX.X 

21  Number  of  layers    XX 


1  An  X  represents  a  space  provided  for  the  variable- 
number  code.  A  number  indicates  a  specific  descriptive 
adjective  to  be  coded  in  the  space  provided  for  that 
variable.  For  example,  a  4  put  in  the  space  designated 
by  the  X  after  "drilling  method"  would  only  select  wells 
drilled  by  an  auger. 
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KTEST  IS  TOTAL  NUMBER  OF 
HYDRAULIC  COEFFICIENTS  USED 

LIST  IS  TOTAL  NUMBER  OF 
LAYER  VARIABLES  CALLED 

<  - 

VARIABLE  READ  FROM  CARD  TYPE  A 

Figure  4. — Flow  chart  for 


of  computer  storage  limitations.  The  codes  are 
stored  by  lithologic  layer  in  the  computer.  Two 
primary  or  one  primary  and  one  secondary 
constituent  describe  each  lithologic  layer.  These 
data  are  also  reduced  to  a  form  that  can  show 
the  hydraulic  properties  (coefficients)  of  the 
aquifer. 

The  lithologic  type  and  color  are  based  on 
log  descriptions  made  in  the  field.  These  logs 
need  to  be  interpreted  accurately  before  the 
data  can  be  standardized  and  stored  as  layer 
data.  For  analytical  interpretation  of  test-hole 


logs,  samples  are  taken  from  holes  cored  ad- 
jacent to  formerly  drilled  test  sites.  Grain  size 
of  each  core  sample  is  visually  described  and 
analyzed  for  comparison  with  descriptions 
based  on  drill  cuttings.  The  median  grain  size 
is  equivalent  to  the  grain  size  having  the  high- 
est percentage  of  occurrence  in  the  drill-cutting 
sample  and  is  termed  "primary  constitutent." 
A  grain  size  representing  a  smaller  percentage 
of  the  sample  is  labeled  "secondary  constitu- 
ent." With  these  comparisons  and  terms,  quanti- 
tative relationships  can  be  transposed  into 
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K  ANO  S  VALUES  AVERAGED  IN 
SATURATED  ZONE  ONLY 




K  AND  S  VALUES  AVERAGED  FOR 

TOTAL  TEST  HOLE 

/ WRITES  OUT  ALL  LAYER  VARIABLES 


WRITES  OUT  ONLY  LAYER 
VARIABLES  SEARCHED 


retrieval  control  cards,  A-N. 


terms  for  coding  the  visual  description  of  drill 
cuttings. 

Hydraulic  properties  of  the  aquifer  are  de- 
termined by  using  a  relationship  between  grain 
size  and  the  hydraulic  coefficients  (fig.  5). 
The  hydraulic  coefficients  (the  coefficients  of 
permeability  and  storage)  are  based  on  field- 
pump  tests  of  selected  alluvial  material,  lab- 
oratory tests  of  the  core  samples,  and  grain 
size-specific  yield  relationships  (2).  The  hy- 
draulic properties  are  coded  by  dividing  the 
primary  and  secondary  grain  sizes  into  four 


identifiable  ranges  (table  3).  Undifferentiated 
sand  is  included  in  a  fifth  range. 

A  computer  program  (appendix  D)  is  used 
to  assign  a  range  number  (1-5)  to  the  layers 
of  each  well  on  the  basis  of  the  coded  lithologic 
description  shown  in  table  4.  Hydraulic-coef- 
ficient values  representing  these  ranges  can  be 
arbitrarily  selected  from  a  graphical  relation- 
ship similar  to  the  one  shown  in  figure  5.  The 
grouping  technique  makes  them  applicable  to 
all  drill-log  data,  with  only  the  values  repre- 
sented by  each  code  being  subject  to  change 
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Table  2. — Data-processing  code  for  tape  storage,  stratigraphic  layers 


Variable 


Term 


Code  form1 


Variable 


Term 


Code  form1 


1  Layer  number    XX 

2  Layer  thickness    XX 

3  Type  of  material:-  column  1,  num-  XXXXX 

ber  of  primary  constituents  (zero 
if  unknown)  ;  column  2,  code  of 
primary  constituent;  column  3, 
code  of  adjective  for  primary 
constituent  (zero  if  primary 
constituent  has  no  adjective)  ; 
column  4,  numeral  1  if  second 
primary  constituent  or  secondary 
constituent  has  no  adjective 
(zero  if  no  secondary  constitu- 
ent), or  code  of  secondary  con- 
stituent if  it  has  an  adjective,  or 
code  of  second  primary  constitu- 
ent with  adjective;  column  5, 
code  of  second  primary  constitu- 
ent with  no  adjective  or  second- 
ary constituent  with  no  adjective 
(zero  if  no  second  primary  con- 
stituent or  secondary  constituent 
exists),  or  code  of  adjective  for 
secondary  constituent,  or  code  of 
adjective  for  second  primary 
constituent. 

Constituents    X 

Clay   0 

Silt    1 

Sand    2 

Gravel    3 

Sandstone    4 

Carbonate    5 

Shale    6 

Evaporite    7 

Undifferentiated  alluvium  •  ■  ■  8 
Undifferentiated  bedrock  ■  ■  -  ■  9 


Type  of  material — Continued  :- 

Log-description  adjective   

Very  fine   

Fine   

Medium   

Coarse   

Lime   

Dolomite   

Halite   

Fossiliferous   

Color  of  material:  major  color  is 
the  left  digit,  minor  color  the 
right  digit  (a  zero  in  either  place 
indicates  no  color).3 

Analytical  method   

Visual  accumulation  tube   

Sieve   

Driller  inspection  (field)   

Microscopic  examination   

Hydraulic  coefficients4 
Transmissibility  (T)   (gpd/ft)  .... 

Values  from  0  to  8,000   

Values  from  8,000  to  30,000  . . 
Values  from  30,000  to  80,000  ■ 
Values  from  80,000  to  150,000 

Permeability  (K)  (gpd/ft2)   

Values  from  0  to  80   

Values  from  80  to  300   

Values  from  300  to  800   

Values  from  800  to  1,500   

Storage  coefficient  (S)  (dimension- 
less  ratio)   

Values  from  0.0000  to  0.0001  . . 
Values  from  0.0001  to  0.0009  . . 
Values  from  0.0009  to  0.0020  ■  • 
Values  from  0.0020  to  0.0040  ■  ■ 


X 
2 
3 
4 
5 
6 
7 
8 
9 

XX 


X 
1 
2 
3 
4 

X 
1 
2 
3 
4 
X 
1 
2 
3 
4 

X 
1 

2 
3 
4 


1  An  X  represents  a  space  provided  for  the  variable-number  code.  A  number  indicates  a  specific  descriptive  ad- 
jective to  be  coded  in  the  space  provided  for  that  variable. 

2  For  example,  a  log  description  of  a  primary  and  a  secondary  constituent  as  "medium  sand,  with  some  fine  sand" 
results  in  a  code  name  of  "sand,  medium,  with  some  sand,  fine,"  or  a  number  code  of  12423.  Again,  a  log  descrip- 
tion of  two  primary  constituents  as  "medium  sand,  and  fine  sand"  results  in  a  code  name  of  "sand,  medium,  and 
sand,  fine,"  or  a  number  code  of  22423. 

3  For  example,  a  sample  with  a  major  color  of  red  but  no  minor  color  would  be  coded  20,  one  with  no  major 
color  but  a  minor  color  of  red,  02;  a  sample  with  a  major  color  of  red  and  a  minor  color  of  brown  would  be  coded 
25,  one  with  a  major  color  of  brown  and  a  minor  color  of  red,  52.  A  code  of  00  indicates  no  color  at  all  was  re- 
corded. 

4  Code  numbers  1-4  represent  a  range  of  values  for  T,  K,  and  S. 
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based  upon  the  geological  materials  being 
investigated. 

Retrieval  of  lithologic  data. — Data  can  be  re- 
trieved in  many  forms,  but  the  most  useful  are 
cross  sections  and  maps  that  define  aquifer 
boundaries  and  the  horizontal  or  vertical  dis- 
tribution of  the  hydraulic  properties  in  the 
aquifer.  The  schematic  diagram  of  the  retrieval 
system  (fig.  6)  shows  the  assignment  and  use 
of  the  ranges  and  corresponding  estimated 
values  of  the  hydraulic  coefficients  for  presen- 
tation in  different  graphical  forms.  Retrieval 
and  testing  of  test-hole,  layer,  and  water-level 
data  are  essential  to  provide  information  neces- 
sary for  selection  and  plotting  of  specified 
ranges  of  the  hydraulic  coefficients  (fig.  2). 
Retrieval  of  test-hole  data  provides  information 
necessary  for  selection  and  plotting  of  the 
ranges  at  a  specified  depth  or  subdatum  eleva- 
tion. Water-level  data  are  used  to  select  ranges 
characteristic  of  each  layer  within  the  satu- 
rated zone.  Testing  on  the  range  number  of  the 
hydraulic  coefficient  permits  the  selection  of 
specified  ranges  of  the  hydraulic  coefficients. 


Table  3. — Lithologic  type  and  associated  range 
of  the  hydraulic  coeff  icients 


Range1  Lithologic  description  of  log- 

1  Silt;3  silt  with  some  very  fine  sand;  very 

fine  sand  with  some  silt;  silt  and  very 
fine  sand. 

2  Very  fine  sand;  very  fine  sand  with  some 

fine  sand;  fine  sand  with  some  very 
fine  sand;  very  fine  sand  and  fine  sand. 

3  Fine  sand;  fine  sand  with  some  medium 

sand;  medium  sand  with  some  fine 
sand;  fine  and  medium  sand. 

4  Medium  sand;  coarse  sand  and  gravel; 

medium  sand  with  some  coarse  sand  and 
gravel ;  coarse  sand  and  gravel  with 
some  medium  sand;  medium  and  coarse 
sand. 

5  Undifferentiated  material  in  ranges  2,  3, 

and  4. 


1  Zero  is  used  when  a  range  cannot  be  assigned  to  a 
layer. 

2  Combination  of  primary  and  secondary  constituents 
that  are  classified  as  indicated. 

3  Silt  and  clay  are  undifferentiated. 


Table  4. — List  of  codes  representing  four  groups  of  relative  permeability 

Sand 

Clay  and  silt  Very  fine  to  _.  ,.  Medium  to 

(1)  fine  Fme  to  medium  coarge 

(2)  (3)  (4) 


10000 
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12324 

22324 
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22524 
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22311 
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21122 
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12425 
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10013 

11122 

21123 

12400 

13000 

23124 

20017 

11123 

22223 

12413 

13012 

22513 

10017 

20019 

12200 

22322 

13125 

10019 

23010 
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21010 
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12322 

10123  21013 

10122 
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12410 

12510 
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11013 

11017 

11019 

13011 

16000 
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Figure  5. — Relationship  between  hydraulic-coefficient  ranges,  median  grain  size,  and 

the  coefficient  of  permeability. 


One  important  aspect  of  the  retrieval  system 
is  the  use  of  water-level  records  for  selecting 
data  in  the  saturated  zone.  A  subroutine  which 
computes  a  weighted  average  of  all  measure- 
ments within  a  selected  time  period  is  used  to 
determine  the  average  elevation  of  the  ground- 
water surface  in  each  well. 

GRAPHICAL  PRESENTATIONS 

Examples  are  given  of  five  basic  types  of 
graphical  presentations  of  data  retrieved.  One 


cross-section  display  and  four  possible  map- 
view  displays  are  shown.  A  flow  chart  of  the 
steps  followed  in  preparing  graphical  displays 
and  the  list  of  variable  names,  with  descrip- 
tions, used  to  plot  selected  data  are  shown  in 
appendix  E. 

Cross  sections  of  drill-hole  logs. — Drill-hole 
data  are  punched  on  cards  and  may  be  used 
as  input  data  for  plotter  programs.  The  con- 
trol-card setup  required  to  retrieve  data  for  the 
plotting  of  cross  sections  is  shown  in  figure  7. 
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Figure  6. — Data-retrieval  system  using  hydraulic  coefficients. 
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Figure  7. — Control  cards  for  retrieval  of  data  for  plotting  cross  sections. 
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r  igure  8. — Cross  section  near  Alex,  Okla. 
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Figure  9. — Control  cards  for  retrieval  of  data  for  mapping  distributions  of  selected 

hole  data. 
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Figure  10. — Distribution  of  total  depth  of  test  holes. 


For  plotting  options  available,  a  description  of 
each  retrieval  control  card  used  is  found  in 
appendix  A. 

A  plotted  cross  section  of  drill  holes  near 
Alex,  Okla.,  is  shown  in  figure  8.  Appendix  F 
shows  the  data  cards  used  to  plot  cross  sections. 
The  variable  name,  type  of  plot,  and  column 
in  which  each  variable  is  coded  are  shown  for 
each  card.  These  data  may  then  be  used  on  the 
plotting  equipment  available  to  the  user. 


Maps  of  distributed  drill-hole  data. — The 
total  depths  of  the  test  holes  were  plotted  as  a 
distributive  map.  The  control  cards  for  the 
plotter  program  are  listed  in  figure  9,  and  a 
map  showing  the  distribution  of  the  total  depth 
of  each  test-hole  log  appears  in  figure  10.  The 
map  shows  the  depth  to  which  lithologic  data 
are  available,  since  some  test  holes  were  not 
drilled  to  bedrock.  Appendix  F  contains  the 
data  cards  used  to  plot  selected  data  maps,  with 


13 


1020.  5. 
B  C 


I  2  6 
ABC 


~50i     502  503 
ABC 


A 


CARD  TYPE  J 


CARD  TYPE  G 


CARD  TYPE  F 


515  516 
0  P 

CARD  TYPE  C 


CARD  TYPE  B 


~5 — r 

D  E 


CARD  TYPE  A 


Figure  11. — Control  cards  for  retrieval  of  hydraulic-coefficient  ranges  at  a  specified 

subdatum. 
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Figure  12. — Map  of  hydraulic  ranges  for  a  specified  subdatum  elevation. 


the  type  of  plot,  variable  name,  and  column 
containing  the  variable  code  indicated.  Any  of 
the  21  test-hole  variables  may  be  retrieved  for 
plotting  as  a  distributive  map.  Appendix  E 
gives  a  list  of  variable  names  and  their  descrip- 
tions for  all  example  plots. 

Maps  of  data  from  specified  elevations. — 
Because  hydrogeologic  data  concern  informa- 
tion that  must  be  gathered  at  some  depth  below 
the  surface  of  the  ground,  a  convenient  method 
of  looking  at  various  levels,  or  subdatum  eleva- 


tions, was  developed  by  the  authors.  Digital 
computer  techniques  are  used  to  develop  a 
planar  surface  of  lithologic  data  for  a  specified 
elevation  above  mean  sea  level.  Figure  11  shows 
the  control  cards  for  the  subdatum  searching. 
A  map  of  the  test  holes  that  have  lithology  with 
hydraulic  coefficients  ranging  from  1  to  5  at 
the  specified  subdatum  is  presented  in  figure 
12.  Appendix  F  contains  the  data  cards  needed 
to  plot  the  map  shown  in  figure  12. 
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Figure  13. — Control  cards  for  retrieval  of  weighted  averages  of  the  hydraulic  co- 
efficients of  all  layers  in  each  hole. 


Maps  of  average  coefficients  of  permeability 
and  storage. — When  the  ground-water  flow 
through  a  specified  cross  section  or  the  amount 
of  ground  water  stored  in  a  region  of  a  water- 
shed is  needed,  it  is  sometimes  desirable  to 
average  the  coefficients  of  permeability  and 
storage  of  the  layers  that  the  test  hole  pene- 
trates. Figure  13  shows  the  control  cards  for 
averaging  these  hydraulic  parameters. 

Those  averaged  values  of  K  (permeability), 
when  displayed  in  map  view  (as  shown  in  fig- 


ure 14),  provide  a  distribution  of  test  holes 
that  are  potentially  high  or  low  in  ground-water 
yield.  A  distributive  map  of  S  (storage  co- 
efficient) may  be  produced  using  the  system. 
It  is  possible  to  determine  suitable  well  field 
locations  and  possible  artificial  recharge  areas 
using  areal  maps  (fig.  1)  developed  from  the 
distributed  test-hole  data. 

The  retrieved  data  cards  used  to  plot  the 
map  shown  in  figure  14  are  presented  in 
appendix  F. 
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Figure  15. — Control  cards  for  retrieval  of  an  isopachous  map  of  lithology  having 

permeability  ranges  1-5. 
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Figure  16. — Isopachous  map  of  hydraulic  properties,  ranges  1-5. 


Isopachous  maps  of  specific  hydraulic  prop- 
erties.— An  isopachous  map  of  all  the  lithologies 
within  an  alluvial  system  having  the  same 
hydraulic  properties  was  also  plotted.  Such  a 
map  is  used  by  planners  when  well  fields  are 
being  installed  or  augmented.  The  amount  of 
ground  water  available  in  an  area  may  be  de- 
termined from  such  a  map  (3). 


A  control-card  setup  for  retrieving  data  for 
an  isopachous  map  of  all  subsurface  materials 
having  a  hydraulic  code  1-5  is  presented  in 
figure  15.  The  isopachous  map  plotted  is  shown 
in  figure  16. 

Appendix  F  shows  the  retrieved  data  cards 
used  to  plot  the  map  shown  in  figure  16. 
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APPENDIX  A -CONTROL-CARD  SETUPS,  CARD  SPECIFICATIONS, 
AND  LIST  OF  VARIABLES  FOR  OPTIONS  A  N 

Control-Card  Setups 
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Figure  A-l. — Retrieval  option  A :  special  drill-hole  data. 
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Figure  A-2. — Retrieval  option  B :  type  of  material. 
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Figure  A-3. — Retrieval  option  C:  color  of  material. 
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Figure  A-4. — Retrieval  option  D :  data  from  specified  subdatum  planes. 
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Figure  A-5. — Retrieval  option  E :  data  from  selected  heights  above  bedrock. 
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Figure  A-6. — Retrieval  option  F :  data  for  logging  a  test  hole. 
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Figure  A— 7. — Retrieval  option  G :  data  from  saturated  zone. 
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Figure  A-8. — Retrieval  option  H :    data  from  specified  depths  below  ground  surface. 
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Figure  A-9. — Retrieval  option  I :  average  hydraulic  coefficients  for  the  saturated  zone. 
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Figure  A-10. — Retrieval  option  J:  selection  of  all  test-hole  and  layer  data. 


Card  Specifications 


CARD  TYPE  A 


-^variable  | 
optToFI 

NWELL 
1-5 

LVAR 
6-10 

LVARl 
11-15 

NUM 
16-20 

NVAR 
21-25 

ICOND 
26-30 

I  COL 
31-35 

I  DAT 
36-40 

la 

I  LOG 
46-50 

ILVL 
51-55 

IGSUB 
56-60 

KTESK 
61-65 

KTESS 
66-70 

KSAT 
71-75 

A 

X 

X 

X 

X 

B 

X 

X 

X 

C 

X 

X 

X 

X 

D 

X 

X 

X 

X 

E 

X 

X 

X 

X 

F 

X 

X 

X 

X 

G 

X 

X 

X 

X 

H 

X 

X 

X 

X 

I 

X 

X 

X 

X 

X 

X 

J 

X 

X 

X 

CARD  TYPE  B 


VARIABLE 

=Z> 
 1 

i— 1 

Csl 

hi 

^  CD 

dz  i_n 

^  r-i 

HI  O 

^  Csl 

Ln 

hi  Ln 

^  Csl 

LD 

HI  CD 

hn 

OFTioir^- — 

CD 

zr)  i_n 

Q_  1 

zd  i— i 
o_  1 

ZD  1 
Q_  r— i 

ZD  1 
Q_  LD 

ZD  I 
O-  i— 1 

ZD  1 
Q_  LD 

A 

B 

X 

X 

C 

X 

X 

D 

X 

X 

E 

X 

X 

F 

X 

X 

G 

X 

X 

H 

X 

I 

X 

X 

J 

X 

CARD  TYPF  C 

r— 1 

Csl 

Ln 

LD 

OO 

Ol 

CD 
«— 1 

■ — i 

Csl 

m 

■CT 
i— 1 

Ln 
i— i 

LD 

.—1 

^VARIABLE 

£1 
ZD 

 1 

-J 

 1  CD 

 I  Ln 

 1  i— H 

 1  CD 

 1  CNI 

— i  Ln 

—1  Csl 

 1  CD 

_J  rn 

 iLn 

— ir-n 

 1  CD 

 I  -=d- 

 i  Ln 

— i  cr 

 1  CD 

 I  Ln 

 l  Ln 

 i  i_n 

 1  CD 

 1  LD 

 i  Ln 

_J  LD 
UJ  1 

— 1  CD 

— I  r-s, 

UJ  1 

 i  Ln 

_j  ix 

UJ  1 

 1  O 

 1  OO 

LU  1 

OPTlBfT^. 

CD 

UJ  LA 
3  I 

LU  i— 1 
3  1 
i— i  LD 

UJ  1 
3:  < — I 

 1  r— | 

LU  1 

3  LO 

LU  f 
3:  i— 1 

—  Csl 

LU  l 

3  LD 

i — i  CXI 

LU  1 
3:-— t 
• — 1  hO 

LU  1 

3:  LD 
' — '  N~\ 

LU  1 
3  <— I 

■ — i 

LU  1 

3:  LD 

LU  1 
3:  i — | 
— >  Ln 

LU  1 
3  LD 

~  Ln 

3  i — 1 
— .  LD 

3  LD 
— .  LD 

3  i — 1 
— ■  rx 

3  LD 
— •  IX 

A 

B 

C 

CARD 

CONTAINS  TEST -HOLE  NUMBERS  TO  BE  SEARCHED. 

D 

E 

MORE  THAN 

ONE  TYPE  C 

CARD  MAY  BE 

USED 

FOR  ANY 

F 

G 

OPTION.    MAXIMUM  OF  1,000 

TEST 

HOLES  MAY  BE 

H 

I 

SEARCHED, 

J 



Figure  A-ll. — Description  of  card  types  A,  B,  and  C. 
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Figure  A-12. — Description  of  card  types  D,  E,  and  F. 
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Figure  A-13. — Description  of  card  types  G,  H,  and  I. 
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Figure  A-14. — Description  of  card  types  J,  K,  and  L. 
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Figure  A-15.— Description  of  card  types  M  and  N. 
List  of  Variables 

Variable 

name  Description  of  variable 

Type  A 

Number  of  test  holes  to  be  searched  for  data 
retrieval. 

Number  of  lithologic-layer  variables  to  be 
retrieved. 

Number  of  selected  values  of  lithologic-layer 

variables  to  be  retrieved. 
Number  of  additional  lithologic-layer  data 

to  be  retrieved. 
Number  of  descriptive  well  data  to  be  re- 
trieved. 

Number  of  lithologic  materials  to  be  re- 
trieved. 

Number  of  lithologic  colors  to  be  retrieved. 
If  a  1  is  punched  in  column  40,  data  are 
retrieved  from  subdatum  planes  (see  card 
type  J).   If  column   40   is  blank,  this 
option  is  bypassed. 


Card 

NWELL  A-J  15 

LVAR  A  15 

NUM  A  15 

LVAR1  C-J  15 

NVAR  A-J  15 

ICOND  B  15 

ICOL  C  15 

IDAT  D  15 


Variable 
name 


INC 


ILOG 


ILVL 


IGSUB 


KTESK 


KTESS 


KSAT 


Option 


E 


Format 


Description  of  variable 


G 


H 


Card  Type  A — Continued 

15  An  integer  punched  in  columns  41-45  is  the 
number  of  increments  to  be  searched 
above  bedrock  (see  card  type  K) .  If 
columns  41-45  are  blank,  this  option  is 
bypassed. 

15  A  1  in  column  50  indicates  a  total  test-hole- 
log  retrieval.  If  column  50  is  blank,  only 
the  selected  layer  variable  is  retrieved. 

15  A  1  in  column  55  calls  water  level  for 
defining  saturated  zone  (see  card  type 
N).  If  column  55  is  blank,  this  option  is 
bypassed. 

15  An  integer  in  columns  56-60  indicates  the 
number  of  feet  below  ground  surface  to 
be  searched  by  increments.  If  columns 
56-60  are  blank,  this  option  is  bypassed. 

15  An  integer  in  columns  61-65  indicates  the 
number  of  permeability  values  assigned 
(see  card  type  L).  If  columns  61-65  are 
blank,  this  option  is  bypassed. 

15  An  integer  in  columns  66-70  indicates  the 
number  of  storage-coefficient  values  as- 
signed (see  card  type  M).  If  columns 
66-70  are  blank,  this  option  is  bypassed. 

15  If  column  75  is  a  1,  K  and  S  values  are 
averaged  for  the  total  well  log.  If  column 
75  is  blank,  K  and  S  values  are  averaged 
only  for  the  saturated  zone. 


IPUNH(l) 
IPUNH(2) 

IPUNH(3) 
IPUNH(4) 
IPUNH(5) 
IPUNH(6) 

Note: 

IWELL(h) 


D,E,G 


Card  Type  B 
15          A  1  in  column  5  causes  retrieved  layer  data 


IVAR(h) 


to  be  punched  on  cards. 
J  15          A  1  in  column  10  causes  the  net  thickness 

of  selected  lithologic  layers  to  be  punched 
on  cards. 

B-J  15         A  1  in  column  15  causes  test-hole  (well) 

variables  to  be  punched  on  cards. 
C,F  15  A  1  in  column  20  causes  lithologic  color  to 

be  punched  on  cards. 
B  15  A  1  in  column  25  causes  lithologic  material 

type  to  be  punched  on  cards. 
I  15  A  1  in  column  30  causes  the  weighted  aver- 

age of  hydraulic  coefficients  to  be 
punched  on  cards. 

A  blank  in  appropriate  columns  will  bypass  punch  output  for  that  option. 
Card  Type  C 

A-J  1615        This  array  is  used  to  select  test-hole  data 

by  numbers.  As  many  as  1,000  drill-hole 
numbers  may  be  retrieved  from  this  array 
selectively. 
Note :    n  indicates  array  size. 

Card  Type  D 

A  815         Data  are  retrieved  from  this  array  for  as 

many  as  8  lithologic-layer  variables. 
These  variables  describe  hydraulic  and 
geologic  characteristics  for  discrete 
layers. 


25 


Variable 
name 


Option  Format 


Description  of  variable 


Card  Type  E 

TVALU(w)  A  6F12.4      Data  are  retrieved  from  this  array  for  as 

many  as  8  test  values  (6  per  card).  One 
or  more  values  may  be  tested  for  each 
variable  indicated  on  card  type  D.  Addi- 
tional card  type  E  is  needed  for  more 
than  6  values. 


Card  Type  F 

IADVR(n)  C-J  815         Data  are  retrieved  from  this  array  for  as 

many  as  8  layer-variable  values  in  addi- 
tion to  the  ones  tested  for  using  card 
types  D  and  E. 


Card  Type  G 

ISTOR(w)  A-J  1615        Data  are  retrieved  from  this  array  for  as 

many  as  21  drill-hole  variables.  An  ad- 
ditional card,  type  G,  is  needed  for  more 
than  16  test-hole  variables. 


Card  Type  H 

MAJTS(n)  B  1015        Data  are  retrieved  from  this  array  for  as 

many  as  10  discrete  types  of  material 
based  upon  size-distribution  analyses. 


Card  Type  I 

MACOL(n)  C  1015        Data  are  retrieved  from  this  array  for  as 

many  as  10  discrete  colors  of  material. 


Card  Type  J 


SUBDA  D  F10.0 

SUBTS  D  F10.0 

S'UBIN  D  F10.0 


Maximum  elevation  to  be  searched  (feet 
above  mean  sea  level)  when  selecting  data 
from  horizontal  subdatum  planes. 

Minimum  elevation  to  be  searched  (feet 
above  mean  sea  level)  when  selecting  data 
from  horizontal  subdatum  planes. 

The  distance  in  feet  to  be  stepped  down  for 
subsequent  subdatum  planes  from 
SUBDA  to  SUBTS. 


Card  Type  K 

AINC(tt)  E  1615        Data  are  retrieved  from  this  array  for  as 

many  as  16  increments  (feet  above  bed- 
rock) for  which  data  are  retrieved. 


Card  Type  L 

AKVAL(w)  I  10F8.0      Data  are  retrieved  from  this  array  for  as 

many  as  10  discrete  permeability  values. 


Card  Type  M 

SVAL(w)  I  10F8.0      Data  are  retrieved  from  this  array  for  as 

many  as  10  discrete  storage-coefficient 
values. 


Variable 
name 


IMO 


IDY 
IYR 
ITIM 

JMO 
JDY 
JYR 
JTIM 


Option 


G,I 


G,I 
G,I 
G,I 

G,I 
G,I 
G,I 
G,I 


Format 


Description  of  variable 


Card  Type  N 


15 


15 
15 
15 

15 
15 
15 
15. 


Month  of  beginning  of  time  period  for  which 
ground-water  levels  are  averaged  when 
data  are  retrieved  from  the  saturated 
zone. 

Day  of  beginning  of  time  period. 
Year  of  beginning  of  time  period. 
Time    (military)    of    beginning    of  time 
period. 

Month  of  ending  of  time  period. 

Day  of  ending  of  time  period. 

Year  of  ending  of  time  period. 

Time  (military)  of  ending  of  time  period. 


APPENDIX  B -FORMATS  USED  FOR  DATA  STORAGE 

INPUT -CARD  FORMAT  USED  TO  STORE  DRILL-HOLE  DATA 
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I  5 
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F  11.3 

F  11.3 
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INPUT -CARD  FORMAT  USED  TO  STORE  LAYER  DATA 


14 


F4.0 


F4.0 


IX,  FS.O 


F  3.0  F  3.0 


F7.0 


F  7.0 


F  10.0 


14 


o 

X 

UJ° 


COLUMNS  53-80 
BLANK 


(I) 


tr 

UJ 

-j 

(2) 


(3) 


(4) 


(5) 


i 

<r 
i- 

(6) 


Z 

(7) 


(8) 


INPUT -CARD  FORMAT  USED  TO  STORE  WATER-LEVEL  DATA 


12 

17 

15 

F  7.2 

17 

15 

F7.2 

17 

15 

F7.2 

WELL  NO.  * 

|  CODE  NO.  »* 

DATE 

MONTH,  DAY,  YEAR 
2  DIGITS  EACH 

MILITARY  TIME 

DEPTH  TO  WATER  ■»** 
FEET  BELOW 
TOP  OF  PIPE 

E 

o 

TIME 

X 

t- 

0. 
UI 

a 

DATE 

u 
2 

F 

DEPTH 

COLUMN  64-60 
BLANK 

*    ALL  I  FORMAT  INPUT  DATA  IS  RIGHT  JUSTIFIED 
»*  CODE  I  •  RECORDER- CHART  DATA  —  CODE  2  ■  TAPE  -  DOWN  DATA 
•»»  NO  ZERO  PUNCHES  NEEDED  IN  BLANKS  ON  CARDS 


Figure  B-l. — Input  card  formats. 


APPENDIX  C.-SOURCE  LISTING  FOR  STORAGE  AND  RETRIEVAL 

PROGRAMS 


CARD 

I 

2 

APPENDIX  C 

3 
^ 

r 

RETRIEVAL  PROGRAM   FOR  H Y DROGEOLOG I C  DATA 

c 

Q 

SUBROUTINES  USED  NAMED    (CODE.    KTESX »   COLOR.   SERCH.   MATRL  »    TAPES*  JULDY» 

*L 

o 

r 

WATLV  > 

7 

f 

DEFINE   FILE   1 ( 1 000 » 1 6 1 . U > I  SEC  )  . 2 ( 2 1  . 1 000 . U , I  SEC  1 ) » 3 < 1 0 » 3 20 . U . 

IM  ) 

g 

DIMENSION   V(21 ) »VV(21 ) >VAR( 8) 

9 

DIMENSION   MAJTS<  10  )  .MACOL  (10)  .WATEH1000)  »AINC(20)  »I  PUNH  (  6  ) 

1  0 

COMMON    IWELL(1000),IVAR18)»TVALU(8)»IADVI8)  .WVAR ( 1000 )  ♦  I  ST OR  < 

21  ) 

X  l 

COMMON  HEAD(2.252).TITLE(2»10S).JWELL(1000)  » I ADVR ( 8 )  »  STR AT  i  20 

.8  ) 

1 2 

COMMON    ISECISEC1.ISEC2.AKVAL(10)  .SVAL (101.KTEST.KSAT.LTEST 

1 3 

c 

HEADERS  READ    IN   FOR   LABELING  VARIALBLES. 

14 

c 

CONTROLS   FOR  SELECTING  WELL  VARIABLES. 

1 5 

FORMAT ! 16  15) 

16 

c 

1 7 

c 

NWELL-NUMBER  OF  WELLS   TO   BE  SEARCHED. 

1  8 

c 

LVAR  -NUMBER  OF   LAYER   VARIABLES   TO  BE  TESTED. 

1 9 

c 

LVAR 1-NUMBER  OF   LAYER    TEST   VALUES.    TVALU    (I).    TO  BE  READ  IN. 

20 

c 

VALUES  OF    (TVALU    (I)    FORMAT    ( 6F 1 2 . 4 » /2 F 1 2 . 4  )  ) 

2  1 

c 

NUM     -NUMBER  OF  ADDITIONAL   LAYER   VARIABLES  NOT   USED  FOR  SEARCHING. 

22 

c 

NVAR   -NUMBER  OF  WELL  VARIABLES 

23 

c 

ICOND-'O'-BYPASS 

24 

c 

'l'-CALLS    IN   SUBROUTINE   FOR    TESTING   TYPE   OF  MATERIAL  IN 

THE  LAYERS 

25 

c 

ICOL-  'O'-BY-PASS 

26 

c 

•l'-CALLS    IN  SUBROUTINE   FOR    TESTING   LAYER   COLORS  * 

27 

c 

IDAT-  'O'-BY-PASS 

28 

c 

'l'-READS    IN   SUBSTRATUM   ELEVATION   TO   BE   TESTED.  (SUBDA) 

29 

c 

INC  -    'O'-BY  PASS 

30 

c 

'N'-READS    IN  NUMBER   OF    INCREMENTS    IN   FEET   ABOVE  BEDROCK 

TO  BE 

31 

c 

SEARCHED. 

32 

c 

ILOG-    'O'-WRITES  OUT  LAYER  VARIABLE  SEARCHED 

33 

c 

'l'-WRITES  OUT   ALL   LAYERS    IN  A  WELL 

34 

99997 

CONTINUE 

35 

c 

ILVL-    '0'    BY  PASS 

36 

C 

'1'    WATER   LEVEL   CALLED   FOR   SEARCHING    IN   SATURATED  ZONE 

37 

C 

IGSUB-'O'    BY  PASS 

38 

C 

'1'    STEPS  DOWN    INCREMENTALLY  BELOW  GROUND  SURFACE 

39 

C 

KTESK-'O'    BY  PASS 

40 

c 

'N'    NUMBER  OF  PERMEABILITY   VALUES   TO  BE  USED 

41 

c 

KTESS-'O'    BY  PASS 

42 

c 

'N'    NUMBER  OF   STORAGE  COEFFICIENT   VALUES   TO  BE  USED 

43 

c 

KSAT-   '0'    BY  PASS 

44 

c 

,     '1'    K  AND  S  VALUES   ARE   AVERAGED   FOR   TOTAL  THICKNESS 

45 

c 

PUNCH  A    '1'    IN    ILVL  AND  A    'C    IN   KSAT   TO  AVERAGE  K  AND  S   VALUES    IN  THE 

46 

c 

SATURATED  ZONE 

47 

c 

FORMAT   FOR  VARIABLES  NWELL   THROUGH  KSAT    IS  1515. 

48 

99998 

CONTINUE 

49 

c 

I PUNH    (N)    -    '0'    NO  PUNCH      '1'  PUNCH 

50 

C 

THE  FOLLOWING  PUNCH  COMMANDS  ARE    IN   SIX   SEPARATE   FIVE  COLUMN 

FIELDS. 

5 1 

c 

THE   NUMBER   FOLLOWING   EACH  COMMAND    IS   THE   FIELD  NUMBER. 

52 

c 

LAYER  DATA  (1) 

r 

NET  THICKNESS  (2) 

54 

Q 

WELL  VARIABLES  (3) 

55 

c 

COLOR  OF   LITHOLOGY  (4) 

56 

c 

MATERIAL   TYPE   OF  LITHOLOGY  (5) 

57 

c 

WGT •   AVE.    PERM.   AND  WGT.    AVE.   STOR.   COEFF.  (6) 

58 

c 

VALUES  OF   SUBDA.   SUBTS.    SUBIN    (FORMAT  (3F10.0)) 

59 

c 

SUBDA-SUBSTRATUM  ELEVATION 

60 

c 

SUBTS-MINIMUM  ELEVATION   TO  BE  SEARCHED 

61 

c 

SUBIN-INCREMENT    IN   FEET   TO  BE   SUBTRACTED  FROM  SUBDATUM 

62 

c 

SUBROUTINE   TAPES  LOADS   TEST   HOLE   AND  LAYER   DATA   FROM   TAPES  TO 

DISC  FILE 

63 

c 

FOR  SEARCHING. 

64 

CALL  TAPES 
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65  DO  1000  J=l»2 

66  C  READ  LAYER  VARIABLE  NAMES   FROM  CONTROL  CARDS 

67  1000  READ ( 5. 106  )  (HEAD  I  J. I ) ♦ 1  =  1  .108  ) 

68  106  F0RMAT(72A1./.36A1) 

69  IM-1 

70  C  WRITES  LAYER  VARIABLE   NAMES  ON  DISC 

71  WR I TE ( 3 ' I M ) HEAD 

72  DO  999  J=1.2 

73  C  READ  TEST  HOLE  VARIABLE  NAMES  FROM  CONTROL  CARDS 

74  999   READ! 5. 107  )  (HEAD (J. I )»  1  =  1  .252  ) 

75  107  F0RMAT(72A1/72A1/72A1/36A1  ) 

76  IM=3 

77  C  WRITES   TEST  HOLE  VARIABLE   NAMES  ON  DISC 

78  WR I TE I  3 ' I M ) HEAD 

79  1012  READ ( 5  »      100 ) NWE LL . LVAR » L VAR 1 t NUM # NVAR • I COND . I  COL » I  DAT • I NC » I  LOG » I L 

80  1VL» IGSUB.KTESK .KTESS.KSAT 

81  C  TEST  ON  NUMBER  OF  TEST  HOLES  TO  BE  SEARCHED 

82  IF(NWELL) 1013.1013.4000 

83  C  READ  TEST   HOLE  NUMBERS   TO  BE   SEARCHED    (1615)    INTO  ARRAY 

84  4000  READ(5.100)  ( IPUNHI  I )  .1  =  1.6) 

8  5  READ ( 5.100) ( I  WELL ( I ) . 1=1 .NWELL) 

86  100  FORMAT(16I5) 

87  C  TEST  ON  NUMBER  OF  LAYER  VARIABLES  TO  BE  TESTED 

88  4001  IF(LVAR)4003»4003.40Q2 

89  C  READ  VARIABLES  TO  BE  TESTED   INTO  ARRAY 

90  4002  READ15. 100)  (  IVARd  )  .1-1. LVAR) 

91  C  TEST  ON  NUMBER  OF  TEST  VALUES  TO  BE   TESTED  FOR 

92  4003    IF.(LVAR1)4005. 4005.4004 

93  C  READ  VALUE  OF  LAYER  VARIABLE   TO  BE   TESTED   INTO  ARRAY 

94  4004  READ( 5 » 101  )  (TVALUt I )  . 1  =  1 .LVAR1 ) 

95  101  FORMAT    ( 6F 12 . 4  . / 2F 12 . 4 ) 

96  C  TEST   ON  NUMBER  OF  ADDITIONAL  LAYER  VARIABLES  NEEDED! 1-7) 

97  4005    IF(NUM)4007, 4007. 4006 

98  C  READ  ADDITIONAL  LAYER  VARIABLES    INTO  ARRAY 

99  4006  READ< 5. 100 >  ( IADVRl I )  • 1=1 »NUM) 

100  C  TEST   ON  NUMBER  OF  ADDITIONAL   TEST   HOLE   VARIABLES  (1-21) 

101  4007    IFINVAR 14009.4009.4008 

102  C  READ  ADDITIONAL   TEST  HOLE  VARIABLES   INTO  ARRAY 

103  4008  READf 5.100) ( ISTORl I ). 1=1 .NVAR ) 

104  C  TEST    IF  SUBDATUM  ELEVATION   IS  TO  BE  SEARCHED 

105  4009   IF( IDAT13021  .  3021  . 3023 

106  C  TEST   IF  A  HEIGHT  ABOVE  BEDROCK   IS  TO  BE  SEARCHED 

107  3021   IF(  INO  3022  . 3022  . 3024 

108  3024  READ! 5  . 11 12  )  ( AINCI I ) . 1=1 . INC ) 

109  1112  FORMAT ( 16F5  •  1 ) 

110  GO  TO  3124 

111  3022  SUBDA=0.0 

112  GO  TO  3124 

113  C  READ  FROM  ONE  CARD. ( 1 ) ELEV  TO  BE  SEARCHED .( 2 ) ELEV  TO  STOP  ON. (3) 

114  C  INCREMENT   TO  STEP  DOWN-FEET 

115  3023  READ15.5001SUBDA.SUBTS.SUBIN 

116  500  FORMAT ( 3F10.0 ) 

117  C  TEST    IF  MAJOR  OR  MINOR   LITHOLOGIC   TYPE  WANTED 

118  3124   IF( ICOND) 5002.5002.5001 

119  C  READ  CODE  OF  LITHOLOGY  WANTED   INTO  ARRAY 

120  5001   READ! 5.100) (MAJTSl I ) .1=1 . ICOND) 

121  C  TEST    IF  MAJOR  OR  MINOR   COLOR  OF  LITHOLOGY  WANTED 

122  5002   IF( ICOL) 5004.5004.5003 

123  5003  RE AD (5. 100) (MACOL( I ) .1=1. ICOL) 

124  5004  L I ST=NUM+LVAR 

125  KTEST=KTESK+KTESS 

126  IF(KTESK) 1122.1122. 1120 

127  1120  READ  (5*1121)  (  AKVAL  (I)»I  =  1.KTESK.) 

128  1121  FORMAT ( 1 0F8 • 0 ) 

129  1122    IF(KTESS) 1124.1124,1123 

130  1123  READt 5.1121) (SVALl I ) .1=1 .KTESS) 

131  1124  CALL   CODE (NWELL. LVAR. LVAR1 .NUM.NVAR . ICOND. ICOL. IDAT. INC. ILOG. IPUNH 

132  1.ILVL.MAJTS.MACOL.LIST.SUBDA.SUBTS.SUBIN. IS TOR. I ADvR »KSAT  .  IGSUB »KT 

133  2ESK. KTESS. AKVAL. SVAL.AINC) 

134  IF( ICOND) 5006.5006.5005 

135  5005  CALL  MATR L ( I COND  .NWELL »MA JTS » I PUNH ( 5 )  ) 
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136 

5006 

IF( ICOL)5080»5080»5007 

137 

5007 

CALL   COLORl  ICQLfNWELL«MAC0L»IPUNH(4)  ) 

138 

5080 

TCltrCIIEIT.l  1 

I F ( I GSUB ) 2  •  2  •  1 

139 

1 

SUBDA  =  1  •  0 

140 

AMIN=9999» 

141 

GO  TO  8002 

142 

2 

t  c  f  c  i  id  r>  a  \  onnA.onnn    q  o.  a  *3 
I  r  lbublJA/oUUU»oUUUtoUUic 

143 

suuu 

t  c  i  t  i  v/i  i  onm  .  Dnn?  .  onn  i 
1  r  \  iLVL ) ouu^  »ouu^  touui 

144 

8001 

^  A  1    1        lil  A  V  1    \  /  /  Ml  f  f"  1    1     .  T  MCI    1     .  i,|  A  T  T  1  \ 

CALL   WA  T L  V  I  NWE  LL  i  I  WE L L  •  W  A  T  E  L  ) 

145 

8002 

IF(LIST)6004»600^t3014 

146 

3014 

I  M=  1 

147 

C 

READS   HEADERS   Or    LAYER   VARIABLES   FROM  DISC 

148 

READ ( 3 1 IM ) HEAD 

149 

DO  1003  K=l»2 

150 

J  =  0 

151 

c 

TtST    FOR   NUMBER   Or    VARIABLES    TO   BE  SELECTED 

152 

IF(LVAR)90j»903i70i 

153 

902 

DO  1002  JC=1»LVAR 

154 

C 

SELECTS  HEADERS  FOR  VARIABLES 

15  5 

ISTRT=(  I VAR ( JC ) *12— 1 1 ) 

156 

IEND=ISTRT+11 

157 

DO  1001  L=ISTRT»IEND 

158 

J  =  J+l 

159 

1001 

TITLE!K»J)=HEAD(K»L) 

160 

1002 

CONT I NUE 

161 

c 

TEST   FOR  ADDITIONAL  LAYER  VARIABLES 

162 

IF (NUM) 1003»1003»903 

163 

903 

DO  2002  JC=ltNUM 

164 

C 

SELECTS  HEADERS   FOR   LAYER  VARIABLES 

165 

ISTRT=t IADVR(JC)*12-11) 

166 

I  END= I STR T+l 1 

167 

DO  2001  L=ISTRT»IEND 

168 

J  =  J+1 

169 

2001 

TITLE<K»J)=HEAD(KiL) 

170 

2002 

CONT I NUE 

171 

1003 

CONT I NUE 

172 

C 

TEST  ON  SUBDATUM  ELEV 

173 

IF(SUBDA) 5016»5016»3012 

174 

5016 

IF(ILVL)3013»3013»5017 

175 

C 

LABELS   TEST  HOLE  N0.»   SUBDATUM.   DEPTH   TO  TOP  OF 

LAYER 

i   ELEVATION  OF  THE 

1 76 

C 

TOP  OF  A   LAYER »   AND  ELEVATION  OF   THE   TOP  OF  SATURATED 

ZONE  • 

177 

5  017 

WRITE(6»109) 

178 

WRITE(6»5018) 

179 

5018 

FORMAT (  '  +  • • 1     WN  SUB       LD  L  ELEV  SATURATED  ZONE  1 

1  ) 

180 

GO  TO  3013 

181 

3012 

WRITE(6tl09> 

182 

C 

LABELS   TEST  HOLE  NO*  t   SUBDATUM*   DEPTH  TO  TOP  OF 

LAYER 

i   ELEV  OF  TOP  OF 

183 

C 

LAYER    IN  WHICH  SUBDATUM    IS  FOUND 

1 84 

109 

FORMAT ( ' 0 1 ) 

lac 

185 

WRITE(6»108) 

186 

108 

FORMAT ( ' + ' » '     WN  SUB       LD  L  ELEV1) 

18  7 

3013 

IEND=J 

188 

DO  1004  K=lt2 

189 

C 

WRITES  HEADERS  FOR  VARIABLES 

190 

1004 

WRITE! 6 »1 15) (TITLE (K»J) »J=1.IEND) 

191 

115 

FORMAT    !  '    1  1 22X*   96A1  ) 

192 

ELV1=0.0 

193 

IC0NT=0 

1 94 

JCONT=0 

195 

C 

SET  BRANCH  «1 

196 

c 

READ  A  TEST  HOLE  NUMBER   TO  BE  SEARCHED 

197 

7001 

SUB=SUBDA 

198 

IC0NT=ICONT+l 

199 

IF!  I NC ) 1 1 1 0  » 1 1 1 0  » 1 1 1 3 

200 

1110 

I  I NC  =  1 

201 

GO  TO  1114 

202 

1113 

I  I NC=  I  NC 

203 

1114 

DO  1111    I  I  =  1  •  I  INC 

204 

DO  10  IC*ltNWELL 

205 

LTEST=0 

206 

IBRN=1 
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207 

I  SEC  =  2 1 

208 

JBRN=2 

209 

5UBDA=SUB 

210 

IF(SUBDA) 5000.5000.3001 

211 

5000 

IF  ( IlVL)3001»3001»502 

212 

502 

JBRN=1 

213 

c 

READS  TOTAL  NUMBER  OF  LAYERS  FOR   EACH  TEST 

HOLE  FROM  DISC 

214 

300 1 

READ   ( 2 ' I  SEC ) WVAR 

215 

I KKa I  WELL ( I C ) 

216 

c 

BRANCHES  TO  DIFFERENT  PARTS  OF  PROGRAM 

217 

GO  TO ( 3002 » 3003  •  53 ) . I BRN 

2  1  ft 

c 

SETS  MAXLA  EQUAL  TO  TOTAL  NO.   OF   LAYERS  FOR 

EACH  TEST 

HOLE 

219 

3002 

MAXLA=WVAR( IKK) 

220 

GO  TO  3004 

221 

c 

READS  GROUND  SURFACE   ELEV  FOR   EACH  TEST  HOLE  FROM  DISC 

222 

3003 

GSELV=WVAR ( IKK ) 

£.  i.  3 

GSURF«GSELV 

224 

c 

SETS  SECTOR  TO  FIND  BEDROCK  ELEV 

225 

ISEC=7 

226 

c 

SET  BRANCH  *  3 

227 

IBRN=3 

228 

GO  TO  3001 

229 

c 

READS  BEDROCK   ELEV  FOR  EACH  TEST  HOLE  FROM 

DISC 

230 

53 

BEDRK*WVAR ( IKK ) 

231 

c 

TEST   IF  SUBDATUM  ABOVE  BEDROCK   IS  DESIRED 

232 

IF(  INO  3006. 3006. 9001 

233 

c 

ESTABLISH  SUBDATUM  ELEV  ABOVE  BEDROCK 

234 

900  1 

SUBDA=BEDRK+AINC (II) 

235 

GO  TO  3006 

236 

c 

TEST    IF  SUBDATUM  ABOVE  BEDROCK  EXISTS 

237 

3004 

IF (SUBDA) 9999*99  99,3005 

238 

9999 

IF( ILVL )9002»9002i3005 

239 

c 

TEST    IF  SUBDATUM  ABOVE  BEDROCK   IS  DESIRED 

240 

9002 

IF ( INC ) 3006.3006  »3005 

241 

c 

SETS  SECTOR  TO  SELECT  GROUND  SURFACE  ELEV. 

242 

3005 

ISEC=6 

243 

C 

SET  BRANCH  =  2 

244 

IBRN=2 

245 

GO  TO  3001 

246 

3006 

IF ( IGSU8 ) 12001 .12001 .2000 

247 

2000 

SUBDA=GSURF 

248 

SUBDA=SUBDA-( IGSUB*ICONT ) 

tt" 

SUBTS=BEDRK 

250 

IF (SUBTS) 2003.2  003.20  05 

251 

2003 

WRITE(6.2004) IWELL( IC) 

252 

2004 

FORMAT! '0' »25( '*') .'NO  BEDROCK  ELEVATION  FOR  WELL  '.14 

»25(  '*•■)  »■//) 

253 

GO  TO  10 

254 

2005 

IF (BEDRK-AMIN ) 2006.12  001 » 12  001 

255 

2006 

AMIN=BEDRK 

256 

12001 

ISEC2= IWELL ( IC ) 

257 

ILAY«1 

258 

READ   ( 1' ISEC2)STRAT 

259 

SUM=0.0 

260 

THIK-O.O 

261 

7005 

IWLL=IWELL( IC) 

262 

c 

TEST    IF   TEST  HOLE  HAS  A  LOG 

263 

7006 

I F ( MAXLA ) 10.10.200 

264 

c 

TEST  HOLES  WITH  MORE   THAN   TWENTY   LAYERS  ARE 

WRITTEN  ON 

PRINTER. 

265 

200 

IF(MAXLA-20)203.203.202 

266 

202 

WRITEI6.104) IWELLI IC) 

267 

104 

FORMAT ( 1    '.'LAYER  NUMBER   GREATER   THAN   20  ON 

TEST  HOLE 

'  .  I  A  ) 

268 

GO  TO  10 

269 

c 

READS   THROUGH  TEST  HOLE  LAYERS 

270 

203 

IF(  ILVD9997. 9997. 5009 

271 

5009 

SUBDA=WATEL (IKK) 

272 

9997 

DO  90  JLAY=1. MAXLA 

273 

c 

SUMS   LAYER  THICKNESS 

274 

SUM=SUM  +  STRAT (  I  LAY  .  2 ) 

275 

c 

TEST    IF  SUBDATUM  ELEV   IS  READ 

276 

5 

CONTINUE 

277 

5008 

IF (SUBDA) 3008.3008.3007 
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278 

3007 

IF(SUBDA-GSELV )4016 .4016 ♦ 10 

279 

4016 

IF(SUBDA-BEDRK )10»4017.4017 

280 

C 

SUBTRACT   LAYER   THICKNESS  OFF  GROUND  SURFACE 

281 

4017 

ELEV=GSELV-STRAT ( I  LAY  ,2 ) 

282 

C 

TEST   NEW  GROUND  SURF   AGAINST   SUBDATUM  ELEV 

283 

IF   ( ELEV-SUBDA)3010»3010.3009 

284 

C 

RESETS  GROUND   SURFACE    TO   TOP  OF  NEXT  LAYER 

285 

3009 

GSELV=ELEV 

286 

C 

INCREASE   LAYER  NO.   BY  1 

287 

I  LAY  = I L AY+ 1 

288 

GO  TO  90 

289 

C 

COMPUTES  DEPTH   TO  TOP   OF   LAYER  CONTAINING  SUBDATUM 

290 

3010 

GO  TO   ( 5010.5012 ) » JBRN 

291 

5010 

JBRN=2 

292 

SUM=GSURF-SUBDA 

293 

STRAT (  ILAY.2)=SUBDA-ELEV 

294 

ELEV=SUBDA 

295 

GSEL V=SUBDA 

296 

SUM=SUM+STRAT (  I  LAY, 2 ) 

297 

IF(KSAT) 3007.3007.7008 

298 

7008 

CALL   K.TESX  I  ILAY.S7.S8  » AK  V  1  »SVl  ) 

299 

5012 

DTOLA=GSURF-GSELV 

300 

C 

TEST    IF   ADDITIONAL   LAYER  VARIABLES 

301 

3020 

IF(NUM)3016. 3016. 3017 

302 

3017 

DO  3015  IJJ=1.NUM 

303 

C 

SETS  K   TO  ADDITIONAL   LAYER  VARIABLE 

304 

K=IADVR< IJJ) 

305 

C 

SETS    IV   TO   THE  NUMBER   OF   THE   ADDITIONAL   LAYER  VARIABLES 

PLUS   THE   NUMBER  OF 

306 

c 

LAYER  VARIABLES   TO  BE  SELECTED 

307 

IV=I JJ+LVAR 

308 

c 

STORES  ADDITIONAL  LAYER   VARIABLES    IN  ARRAY 

309 

3015 

VAR( IV)=STRAT( ILAY.K) 

310 

IF1LVARI3016. 3025. 3016 

311 

3016 

DO  3011  IJ=1.LVAR 

312 

C 

SETS  <  TO  LAYER  VARIABLE   TO  BE  SELECTED 

313 

K=IVAR( IJ) 

314 

C 

STORES  LAYER  VARIABLES  SELECTED   IN  ARRAY 

315 

3011 

VAR( IJ)=STRAT( ILAY.K) 

316 

C 

COMPUTES  DEPTH   TO  TOP   OF    LAYER  WHICH  CONTAINS  SUBDATUM 

317 

3025 

ELV=GSURF-(SUM-STRAT( ILAY.2) ) 

318 

C 

SETS    IDT   TO  DEPTH  TO  LAYER 

319 

I DT  =  DTOLA 

320 

C 

SETS    ISUB    TO  SUBDATUM  ELEVATION 

321 

ISUB=SUBDA 

322 

IF(KSAT)300. 300. 7003 

323 

C 

TEST    IF   LOG  OF  ALL  LAYERS   IS  WANTED 

324 

300 

IF( ILOG17050. 7050. 7003 

325 

C 

WRITE  DATA  OUT   FOR  ALL  LAYERS 

326 

7050 

IF( IGSUB)7051  . 7051  , 7003 

327 

C 

LABELS  SUBDATUM  ELEVATION,    DEPTH   TO   LAYER,    ELEVATION  OF 

TOP  OF 

LAYER 

328 

c 

CONTAINING  SUBDATUM,   AND   SELECTED  VARIABLES. 

329 

7051 

IF (ELV-ELV1 ) 7003 ,10 ,7  003 

330 

7003 

WRITE    (6,113)    IWLL,            I  SUB . I DT . ELV . ( VAR < I  )  .  I  =  1 » L I S T ) 

331 

113 

FORMAT (  1    •  .  13.  I5.I4.F8.2  .8F12.4) 

332 

C 

TEST   FOR  PUNCHED  OUTPUT. 

333 

IF( IPUNHl 1 ) )94.94,1050 

334 

1050 

WRITE(7,1051)IWLL.ISUB.IDT.ELV.( VAR (I), 1  =  1, LIST) 

335 

1051 

FORMAT ( I3.I5.I4.F8.2.5F12.4./.5F12.4) 

336 

94 

IF(KSAT)301»301»93 

337 

301 

IF( ILVL) 10.10.95 

338 

95 

I LAY= I LAY+ 1 

339 

C 

TEST   LAYER   VARIABLE  VALUES  SEARCHED. 

'40 

C 

LVAR1    IS   THE  NUMBER  OF  VALUES  SELECTED   FOR   EACH   LAYER  VARIABLE 

TESTED. 

341 

3008 

IF1LVAR1 )4018. 4018, 904 

342 

904 

DO  6  JJ=1.LVAR 

343 

K=IVAR( JJ) 

344 

DO  6  J=1.LVAR1 

345 

C 

TEST   DATA    IN  LAYER  AGAINST   TEST  VALUE 

346 

I F  t  TVALU ( J ) -STRAT (JLAY.K)  16.4011.6 

347 

6 

CONTINUE 

348 

GO  TO  9 
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349 

C 

TEST   FOR  ADDITIONAL  LAYER  VARIABLES 

350 

4011 

THIK=THIK+STRAT( JLAYf 2 ) 

351 

4018 

IF ( NUM  >  6003  1 6003 ,4010 

352 

4010 

DO  8    U  =  1»NUM 

353 

K= I ADVR ( IJ) 

354 
355 

C 

SUMS  NUMBER  OF  LAYER  VARIABLES  TESTED  AND  ADDITIONAL  LAYER 
IV* I J+LVAR 

VAR I ABLES 

356 

8 

VAR(  I  V ) =STRAT ( JLAY »K ) 

357 

C 

TEST  ON  NUMBER  OF  LAYER  VARIABLES  TO  BE  SEARCHED 

358 

6003 

IF ( LVAR ) 6001 ♦ 600 1 » 6000 

359 

6000 

DO  7  U=ltLVAR 

360 

C 

SETS  K  TO  THE  VARIABLE  TESTED 

361 

K=IVAR(  I  J) 

362 

C 

SETS  VAR(IJ)    TO  THE  VALUE  STORED   IN  LAYER  DATA 

363 

7 

VAR< I J)=STRAT( JLAY.K) 

364 

C 

TEST  ON  VARIABLES  TO  BE   SEARCHED  FOR  OR  ADDITIONAL 

LAYER  VARIARLFS 

365 

6001 

IF(LIST)6004»6004»6002 

366 

6002 

IF(  ILVD6020.6020.9 

367 

6020 

WRITE (6. 103)  IWLL»JLAY> ( VAR (  I  )  »  I  =  1 tL I  ST ) 

368 

103 

FORMAT ( '    ' »10X.2I5.8F12.4) 

369 

C 

TEST   FOR  PUNCHED  OUTPUT 

370 

IF(  IPUNHI 1 )  )9»9«6006 

371 

6006 

WRITE(7»6007) I WLL  » JLAY  » (VAR( U) »IJ=1»LIST) 

372 

6007 

FORMAT(2l4,8F9.2/»8X,8F9.2/»8X»5F9.2) 

373 

C 

RETURN   TO  READ  NEXT  LAYER 

374 

9 

IF(KTEST)90»90»6008 

375 

6008 

IF( ILVL)92»92»91 

376 

91 

KLAY= I LAY-1 

377 

GO  TO  920 

378 

92 

KLAY* JLAY 

379 

C 

SUBROUTINE  KTESX  COMPUTES  AVERAGE  PERMEABILITY  AND 

STORAGE 

COEFFICIENT 

380 

C 

VALUES. 

381 

920 

CALL  KTESX I KLAY»S7,S8 ,AKV1 ,SV1 ) 

382 

90 

CONTINUE 

383 

93 

IF (KTEST ) 6011  ,6011 ,6009 

384 

6009 

AKV1=AKV1/S7 

385 

SV1=SV1/S8 

386 

WRITE (6, 60 10) I WLL  »AKV1 ,SV1 

387 

6010 

FORMAT (  '    ','HOLE  NUMBER*    ',14,/,'    WGT.   AVE. PERM.*  ' 

,F10.0,/ 

,'  WGT, 

388 

1AVE.   STORAGE  COEFFICIENT*  '.F10.3) 

389 

IF( IPUNHI 6) 16011,6011,11100 

390 

11100 

WRITE! 7,11101 )  IWLL.AKVl »SV1 

391 

11101 

FORMAT ( I4,F5.0,F5.3) 

392 

C 

RETURN  TO  READ  NEXT   TEST  HOLE 

393 

6011 

IFILVARl  1 10,10,4019 

394 

4019 

WRITE! 6 ,161 ) I WLL ,TH IK , ( TVALUI  J) , J  =  l  ,LVAR1 ) 

395 

161 

FORMAT ( 1            HOLE  NO.    ',I5»'    NET   THICK    ',F12.0,'    CODES  '»5X 

,6F4.0) 

396 

IF{ IPUNH(2)  )10,10,4020 

397 

C 

PUNCHES  DATA  ON  CARDS 

398 

4020 

WRITE ( 7 , 161 ) I WLL  »TH IK , ( TVALUI J) , J  =  l ,LVAR1  ) 

399 

10 

CONTINUE 

400 

111.1 

CONTINUE 

401 

c 

TEST    IF  SUBDATUM  ELEVATION   IS  SELECTED 

402 

6005 

IF( ILVL)99,99,98 

403 

98 

SUBDA=0.0 

404 

99 

IF (SUBDA) 6  004,6004,7000 

405 

C 

TEST    IF  POINT  ABOVE   BEDROCK   IS  SELECTED 

406 

7000 

IF(INC)900  3,9003,6004 

407 

c 

STEP  DOWN  SUBDATUM  TO  BE  SEARCHED 

408 

9003 

IF(IGSUB)9004»9004,9005 

409 

9004 

SUBDA=SUBDA-SUBIN 

410 

IF (SUBDA-SUBTS 16004,6004, 7002 

411 

7002 

ELV1=ELV 

412 

GO  TO  7001 

413 

C 

TEST    IF  SUBDATUM  WANTED   IS  BELOW  LOWER  LIMIT 

414 

9005 

IF (SUBDA-AMIN 1 6004 , 6004 , 700 1 

415 

C 

RESET   ELEVATION  OF  TOP  OF  LAYER 

416 

C 

TEST    IF  ADDITIONAL   TEST   HOLE  VARIABLES  ARE  SELECTED 

417 

6004 

IF(NVAR) 1011,1011,11 

418 

C 

SET  DISC  SECTOR 

419 

11 

IM=3 
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420 

C 

READ  HEADERS  FOR  ADDITIONAL   TEST  HOLE  VARIABLES 

421 

READ   (3'IM)  HEAD 

422 

JSTRT=1 

423 

JEND=9 

424 

DO   1009  IL=1»3 

425 

c 

TEST    IF  NUMBER    IS  MORE   THAN  NINE 

426 

IF    (NVAR-JEND) 1010,1010.12 

427 

c 

SETS  NUMBER  OF  VARIABLES   TO  BE  WRITTEN 

428 

1010 

JEND=NVAR 

429 

12 

DO   1007  K=1.2 

430 

J  =  0 

431 

C 

READS  NINE   HEADERS   FROM   DISC  STORAGE 

432 

DO  1006  JC= JSTRT » JEND 

433 

ISTRT=( ISTORI JC>*12-11 ) 

434 

IEND=ISTRT+11 

435 

DO  1005  L= ISTRT  » I  END 

436 

J  =  J+1 

437 

1005 

TITLE ( <» J ) =HEAD( K»L> 

438 

1006 

CONTINUE 

439 

1007 

CONTINUE 

440 

IEND=J 

441 

DO  1008  K=l»2 

442 

C 

WRITES  NINE  HEADERS  ON   EACH  PRINTER  LINE 

443 

1008 

WRITE    (6.105)    (TITLE    (K.J)  » J= 1 » I  END ) 

444 

105 

FORMAT  (  '  '.108AD 

445 

C 

TEST    IF   SUBDATUM  ELEVATION   IS  SELECTED 

446 

DO  22  J=1.NWELL 

447 

DO  20  I=JSTRT»JEND 

448 

c 

STORES  SELECTED  VARIABLES  ON  DISC 

449 

ISEC1=IST0R( I ) 

450 

c 

READS  VARIABLES  FROM  DISC 

451 

READ   12'  ISECnWVAR 

452 

c 

SETS   IK   TO  TEST  HOLE  NUMBER 

453 

I K= I  WELL ( J ) 

454 

c 

SETS    IKK   TO  VARIABLE 

455 

IKK= ISTOR (  I  ) 

456 

V(  IKK ) =WVAR ( IK  ) 

457 

20 

CONTINUE 

458 

DO  21  I=JSTRT.JEND 

459 

K=ISTOR ( I  ) 

460 

21 

VV( I )=V(K) 

461 

C 

WRITES  VARIABLES  FOUND  AT  SUBDATUM  ON  PRINTER 

462 

WRITE    (6.5400)    ( VV I  I )  » I =JSTRT . JEND > 

463 

5400 

FORMAT ( 1  '.9F12.3) 

464 

C 

TEST   FOR  PUNCHED  OUTPUT 

465 

IF(  IPUNHO)  )22  .22.  1054 

466 

1054 

WRITE(7.1055)(VV(I) . I = JSTRT . JEND ) 

467 

1055 

FORMAT(6F12.3./.6F12.3./.6F12.3./.3F12.3) 

468 

C 

PUNCHES  VARIABLES  FOUND  AT   SUBDATUM .   SIX   PER  CARD 

469 

22 

CONT INUE 

470 

C 

TEST    IF   ALL  VARIABLES  SELECTED  HAVE   BEEN  WRITTEN 

471 

IF   ( JEND-NVAR ) 23  .1011 .1011 

472 

c 

RESET   TO  NEXT   NINE   HEADINGS  AND  VARIABLES 

473 

23 

JSTRT=JEND+1 

474 

JEND= JEND+9 

475 

1009 

CONTINUE 

476 

1011 

CONTINUE 

477 

GO  TO  1012 

478 

1013 

STOP 

479 

END 

480 

C 

481 

C 

482 

C 

SUBROUTINE  CODE  DISPLAYS  A  LIST  OF  OPTIONS  SELECTED  BY  THE 

USER. 

483 

c 

RETRIEVAL   CONTROL  CARD  CODING    IS  DISPLAYED  ON  PRINTER 

FOR 

EACH  JOB 

484 

SUBROUT INE  CODE ( NWELL .LVAR . LVAR1 »NUM .NVAR . 1 COND . 

ICOL. 

IDAT  . 

INC. ILOG 

4P5 

1  .  I PUNH  » ILVL  »MAJTS»MACOL » I L I  ST • SUB DA . SUBTS • SUB  I N » 

ISTOR 

.  IADVK.KSAT  » I 

486 

2GSUB.KTESK.KTESS . AK VAL . SV AL . A  I NC ) 

487 

DIMENSION   AKVAL ( 10) »SVAL( 10) iAINC(20) 

488 

DIMENSION  MAJTS( 10) .MACCL ( 10 ) * ISTOR ( 21 ) • I ADVR ( 8 ) 

.  I  PUNH I  6 ) 

489 

COMMON   IWELL(IOOO). IVAR18 )  .TVALU(8) 

490 

WRITEI6.120) 
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491  120   FORMAT!  'l'»'    IF   TEST  CODES  ARE   EQUAL   TO   1,    THEN   TEST    IS  USED.    IF  E 

492  10UAL   TO  0»    THEN   TEST    IS  BY-PASSED  '  »  /  /  ) 

493  WRITE(6il00)  NWELL 

494  100   FORMAT  I  '    '  »  'NUMBER  OF   TEST   HOLES   TO  BE   SEARCHED=  '.13) 

495  WRITE (6.  133)  ( I  WELL (  I )  ♦ 1  =  1  .NWELL) 

496  133  FORMAT (  •    '.'TEST  HOLES   THAT  ARE  BEING   TESTED  ARE'  »/. 41(2415./)  ) 

497  WRITE<6.101 ) LVAR 

498  101   FORMAT (' 0 '.' NUMBER  OF   LAYER  VARIABLES  TO  BE   TES TED=  ',12) 

499  IF  (LVAR12.2.1 

500  1   WRITE ( 6 » 102 )  (  I VARl I )  . 1  =  1  .LVAR  ) 

501  102   FORMAT ( '    '.'THE  VARIABLES   THAT  ARE   TO  BE   TESTED  ARE'/. 815) 

502  2  WRITE16.1031LVAR1 

503  103  FORMAT (' 0 '»' NUMBER  OF   TEST  VALUES=  '.12) 

504  IFILVARD4.4.3 

505  3  WRITE16.104) (TVALU! I ) .I=1.LVAR1) 

506  104  FORMAT (  1    '.'TEST   VALUES  OF   EACH  LAYER  ARE '  •  / » 8F 10 . 4 ) 

507  4  WRITE(6.105)NUM 

508  105   FORMAT! '0' . 'NUMBER  OF  ADDITIONAL   LAYER  VAR I ABLES=  '.12) 

509  IF(NUM)6.6.5 

510  5  WRITEI6.106)  ( IADVR! I )  »I  =  1,NUM> 

511  106   FORMAT ( '    '.'ADDITIONAL   LAYER   VARIABLES  ARE'»/»8I5) 

512  6  WRITEI6.116) I L I  ST 

513  116   FORMAT !  ' 0 '»' TOTAL  NUMBER   OF   LAYER   VARIABLES   THAT  ARE   CAL  LED=  '.12) 

514  WRITE(6.107)NVAR 

515  107   FORMAT (' 0 '»' NUMBER  OF   TEST   HOLE  VAR I ABLES=  '.121 

516  IF(NVAR)8.8»7 

517  7  WRITE(6»108)  ( 1ST OR (  I )  .1  =  1.  M VAR) 

518  108   FORMAT (  '    '.'THE   TEST  HOLE  VARIABLES  ARE'./. 2115) 

519  8  WRITE16.109) IDAT 

520  109   FORMAT (' 0 '.' SU3DATUM  ELEVATION  CODE  =  '.12) 

521  IF( IDAT) 10.10.9 

522  9  WRITE(6.110)SUBDA.SUBTS.SUBIN 

523  110   FORMAT!  '    '.'ELEVATION   TO   BE   SEARCHED=    '.F10.0./.'    ELEVATION   TO  TER 

524  1M I  NATE   SEARCHING^    '.F10.0./.'    STEP-DOWN    INCREMENT*  '.F10.0) 

525  10   IF!  INC)  132  .132  .130 

526  130  WRITE ( 6.131 ) INC» (AINC ( I ). 1=1 » INC ) 

527  131   FORMAT (' ONUMBER  OF   HEIGHT    INCREMENTS  ABOVE   BEDROCK  TO  BE  SEARCHED 
526  1=    '.12./.'    THE    INCREMENT   VALUES  ARE '  .  /  ♦  20F 1 0 . 2 ) 

529  132  WRITEI6.112) ICOND 

530  112   FORMAT (' 0 '»' TEST    IF  MAJOR   OR   MINOR   LITHOLOGY   TYPE  WANTED=  '.12) 

531  IF! ICOND) 12.12»11 

532  11  WRITE(6.113) (MAJTS! I ) .1=1, ICOND) 

533  113   FORMAT!'    '.'LITHOLOGY   TYPE  CODES  ARE'./. 815) 

534  12  WRITE(6.114) ICOL 

535  114  FORMAT (' 0 '»' TEST    IF  MAJOR   OR  MINOR  COLOR  OF   LITHOLOGY  WANTED* ' • 1 3 ) 

536  IF! ICOL) 14,14,13 

537  13  WRITE(6.115) (MACOL! I ) .1=1. ICOL) 

538  115   FORMAT!'    '.'LITHOLOGY   COLOR  CODES  ARE'./, 815) 

539  14  WRITE(6.117) ILVL 

540  117   FORMAT!  '0'. 'TEST    IF  WATER LEVEL  OF   SATURATED  ZONE    IS   TO  BE  USED=  ', 

541  112) 

542  WRI TE ( 6 . 12 1 ) KSAT 

543  121   FORMAT ( 1 OTEST    IF   SATURATED  ZONE    IS   TO  BE   LQGGED=  '.12) 

544  WRITE ( 6 » 118 ) ILOG 

545  118   FORMAT (' 0 '»' TEST    IF  ALL   LAYERS  ARE   TO  BE   LOGGED=  '.12) 

546  IF! IGSUB) 124.124.122 

547  122  WRITE! 6.123 ) IGSUB 

548  123   FORMAT (' OSTEP   DOWN    INCREMENT   FROM  GROUND  SURFACE  ELEVAT I ON=  '.12. 

549  1  '    FEET '  ) 

550  124    I F ( KTESK ) 129.129.125 

551  125  WRITE!6»126)KTESK 

552  126   FORMAT !' ONUMBER  OF  K-PERMEAB I L I T Y  AND  S-STORAGE  COEF.   VALUES=  '.12 

553  1  ) 

554  WRITE(6»127) (AKVAL! I).I=1.KTESK) 

555  127   FORMAT (  '    '.'THE   K-PERMEAB I L I TY  VALUES  ARE '  . / » 10F8 . 0 ) 

556  WRITE(6»128)(SVAL(I).I=1.KTESS) 

557  128   FORMAT!'    '.'THE   S-STORAGE   COEF.   VALUES  ARE '  » / .  1 0F8 . 3 ) 

558  129  WRITE16.119) ( IPUNH! I ) .1=1 .6) 

559  119   FORMAT !' 0 '.' TEST   FOR   PUNCHING   DATA './/. 5X  .' LAYER  DATA=  '»I2»/»5X.' 

560  1NET   THICKNESS*    '  .  I  2 . / ♦ 5X • TEST  HOLE  DATA=    '  .12./ »5X»' COLOR  OF   L  I  THOLCGY  = 

561  1    '  .12./. 5X. 'MATERIAL  TYPE  OF  L I THOLOGY=    '.I2./.5X.'    WGT.   AVE.  PERM 


36 


562 

2.   AND  WGT.   AVE*   STOR.   COEF.=    '  »I2) 

563 

RETURN 

564 

END 

565 

C 

566 

C 

567 

C 

SUBROUTINE   KTESX  COMPUTES  AVERAGE   K  AND  S  VALUES  WEIGHTED  ON  LAYER 

568 

c 

THICKNESS 

569 

SUBROUTINE  KTESX(ILAYiS7iS8»AKV1»SV1) 

570 

COMMON    IWELL(IOOO)  *  1 VAR ( 8  )  »TVALU(8)  t I ADV ( 8 )  »WVARI10O0)  » 1 STOR  t  21 ) 

571 

COMMON  HEAD(2»252) tTITLE(2»108) • JWELL ( 1000)  » I ADVR ( 8 )  iSTRAT120»6) 

572 

COMMON   I  SEC  •  1 5EC1 • I SEC2  tAKVAL ( 10 ) »SVAL ( 10 ) » KTEST tKSAT ♦ LTEST 

573 

I F ( LTEST ) 10  » 10  •  1 5 

574 

10 

S7=0 .0 

575 

S8  =  0  .0 

576 

AKV1 =0 • 0 

577 

SV1=0«0 

578 

LTEST=LTEST+1 

579 

15 

J7  =  STRAT (  I  LAY  »7 ) 

580 

J8  =  STRAT (  I  LAY  »  6 ) 

581 

IF(J7)21»21»20 

582 

20 

S7=S7+STRAT( ILAY»2) 

583 

AKV1=AKV1+(AKVAL(J7)*STRAT( I  LAY • 2 ) ) 

584 

2  1 

IF( J8 ) 26«26»25 

585 

25 

S8=S8+STRAT(  ILAY»2) 

586 

SV1  =  SV1+(SVAL(J8)*STRAT(  ILAY»2)  ) 

587 

26 

RETURN 

588 

END 

589 

c 

590 

c 

591 

SUBROUTINE   COLOR( ICOLRiNWELL.MACOL.IPUNH) 

592 

c 

SUBROUTINE   TO  READ  AND  TEST  FOR  MAJOR  OR  MINOR  COLOR  OF   L I THOLOGY 

593 

DIMENSION  MACOL(IO) 

594 

COMMON    lWELL(100Q)tIVAR(8)  iTVALU(8)»IADV(8)  iWVAR(1000)»ISTOR(21) 

595 

COMMON  HEAD(2»252) »TITLE(2»108) iJWELLl 1000) »IADVR(8) «STRAT(2Q»8) 

596 

COMMON    I  SEC » I  SEC  1 » I SEC2 

EOT 

597 

100 

FUR M ATI"       i ' HOLE   NO.    »2X» 'LAYER   NU«^2X» 'MAJOR   COLOR' »2X»'MINOR  CO 

598 

1LOR 1  »4X < ' ELEV    1  » 3X • 1 DEP T H 1  • 3X i • TH I CKNE SS 1  » 5 X  •  ' NET   THICK(MAJOR)  NE 

599 

2T   TH I CK ( M I  NOR  )  1  ) 

600 

101 

FORMAT  I  1    1  »2X»I3»8X»I2»10X»I2»20X»F6.1f2XF6.1»6X»F6.2) 

601 

102 

FORMAT! 1    1 »2X»I3»8X»I2»23X»I2t7X«F6»l»2X»F6«l»6X»F6»2) 

602 

103 

FORMAT! '+' »80XtF6»2) 

603 

104 

FORMAT! ■ + ' »99X»r6«2) 

604 

c 

WRITE   HEADERS   FOR  COLORS   FOUND   IN  LITHOLOGY 

605 

WR I TE ( 6  » 1 00 ) 

606 

c 

ICOLR   IS   THE  NUMBER  OF  COLORS  SELECTED 

607 

DO  6   K=l . ICOLR 

608 

c 

MACOL(K)    IS  THE  NUMERICAL  CODE  FOR  EACH  COLOR  SELECTED 

609 

JTEST=MACOL ( K ) 

610 

c 

NWELL    IS   THE  NUMBER  OF   TEST   HOLES   TO  BE  SEARCHED 

611 

DO  5  IC=1»NWFLL 

612 

c 

CALLS  SUBROUTINE  NAMED    ' SERCH ' 

613 

CALL   SERCH ( NLAY.GSELV* IC ) 

614 

c 

SETS  SECTOR   TO  TEST  HOLE   NO.   BEING  SEARCHED 

615 

I SEC2= I  WELL ( I C ) 

616 

c 

READS  FROM  DISC   LAYER  DATA 

617 

READ    ( 1 1  I SEC2 ) STRAT 

618 

SUM=  0 . 0 

619 

SUM2  =  0 . 0 

620 

SUM3=0.0 

621 

J=l 

622 

c 

SETS    IWLL    TO  TEST  HOLE   NO.    BEING  SEARCHED 

623 

52 

IWLL=IWELL( IC) 

624 

c 

READS  LAYER  DATA  FOR   ALL   LAYERS    IN   EACH   TEST  HOLE 

625 

53 

DO  4   JJ=1 »NLAY 

626 

c 

SETS  TYPE   TO  COLOR  OF  LAYER 

62  7 

50 

TYPE  =  STRAT ( J«4  ) 

628 

c 

SETS   THICK.   TO   THICKNESS  OF   LAYER   BEING  SEARCHED 

629 

TH  ICK  =  STRA  T ( J  »  2 ) 

630 

c 

SUMS   LAYER   THICKNESS  OF  ALL  LAYERS  AS  SEARCHED 

631 

SUM=SUM+STRAT( J»2) 

632 

c 

KTEST    IS  A  THREE  DIGIT  NUMBER    (COLOR  CODE)   OF  LAYER 
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633  ktest=  i f  i x ( type  ) / 10 

634  c  test  if  color  in  layer  is  selected  color 

635  if  (ktest-jtest)  2.  1»  2 

636  c  sum1  is  the  thickness  to  top  of  layer  where  color  is  found 

637  1  sum1*sum-strat ( j. 2) 

638  sum2=sum2+thick 

639  c  computes  elevation  of  top  of  layer  where  color  is  found 

640  elev=gselv-sum1 

641  c  icol  is  two  number  code 

642  icol  =  type 

643  c  labels  hole  no.,  layer  no.*  color  code.  elev  of  top  of  layer »  depth  to 

644  c  layer  ♦  layer  thickness 

645  write(6»101) iwll.  j j . i  col , ele v t sum  1 • th i ck 

646  c  test  for  punched  output 

647  IF< IPUNH14.4.10 

648  C  PUNCHES  HOLE  NO.  •   LAYER   NO..    COLOR  CODE.   ELEV.   OF   TOP   OF   LAYER.  DEPTH 

649  C  TO  LAYER.   LAYER  THICKNESS  ON  ONE  CARD 

650  10  WRITE(7»11 ) IWLL.JJ. I  COL. ELEV. SUM  1  .THICK 

651  11   FORMAT (  1    COLORl    1 .3  1 5 ♦ F6 . 1 » 2X  » F6 . 1 . 6X . F6 . 2 ) 

652  GO  TO  4 

653  2  KTEST=TYPE-FLOAT(KTEST)*10. 

654  C  TEST   IF  COLOR  ASKED  FOR   IS  MINOR  COLOR 

655  IFIKTEST-JTEST ) 4*3*4 

656  C  CALCULATES  DEPTH  TO  TOP  OF  LAYER 

657  3  SUM1=SUM-STRAT ( J. 2) 

658  SUM3«SUM3+THICK 

659  C  CALCULATES  ELEV.   TOP  OF  LAYER 

660  ELEV=GSELV-SUM1 

661  C  ICOL   IS  A   TWO  NUMBER  CODE 

662  ICOL=TYPE 

663  C  LABELS  HOLE  NO..   LAYER  NO..   COLOR  CODE.    ELEV  OF   TOP  OF  LAYER.  DEPTH 

664  C  TO  LAYER.   LAYER  THICKNESS 

665  WRITEI6. 102) IWLL.  J J . I  COL » ELE V . SUM  1 . TH I CK 

666  IF( IPUNHI4.4.12 

667  C  PUNCHES  HOLE  NO..   LAYER   NO..    COLOR  CODE.    ELEV  UF   TOP  OF   LAYER.  DEPTH 

668  C  TO  LAYER.    LAYER   THICKNESS  ON  ONE  CARD 

669  12  WRITE17. 13) IWLL.JJ. ICOL. E LEV . SUM  1  ♦  TH I CK 

670  13  FORMAT (    'COLOR2  '.3I5.F6.1.2X.F6.1.6X.F6.2) 

671  4  J=J+1 

672  I F ( SUM2  )  55.55.54 

673  54  WRITE(6.103)SUM2 

674  55  IF(SUM3)5.5.56 

675  56  WRITE16.1041SUM3 

676  5  CONTINUE 

677  6  CONTINUE 

678  RETURN 

679  END 

680  C 

681  C 

682  SUBROUTINE  SERCH INLAY .GSELV. IC) 

683  C  SUBROUTINE   READS  FROM   DISC    (1)    TOTAL  NO.   OF   LAYERS.    (2)  GROUND 

684  C  SURFACE  ELEVATION  FOR   SELECTED  TEST  HOLE 

68  5  COMMON    IWELL(IOOO), IVAR ( 8 ) »TVALU( 8  )  » IADV1 8 ) .WVAR (1000>»ISTOR121) 

686  COMMON  HEAD ( 2 .252 ) »TITLE(2»108).JWELL(1000).IADVR(8)»STRAT(2Q.S) 

687  COMMON   I  SEC  .  I  SEC  1 . 1 SE C2 

688  C  SET  SECTOR  NUMBER  TO  READ  TOTAL  NO.   OF   LAYERS  FROM  DISC 

689  IBRN=0 

690  ISEC=21 

691  3001   READ ( 2 1  I  SEC ) WVAR 

692  I KKK  = I  WELL (  IC) 

693  IF( IBRN ) 10.10.3002 

694  10  NLAY=WVAR ( IKKK ) 

695  C  SET  SECTOR  NUMBER  TO  READ  GROUND  SURFACE  ELEV.   FROM  DISC 

696  ISEC«6 

697  IBRN=1 

698  GO  TO  3001 

699  C  SETS  GSELV  TO  GROUND  SURFACE  ELEVATION 

700  3002  GSELV=WVAR ( IKKK) 

701  RETURN 

702  END 

703  C 
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704 

C 

705 

SUBROUTINE   MATRL    t I COND . NWELL » MA J TS » I PUNH ) 

706 

C 

SUBROUTINE  SELECTS  MAJOR   AND  MINOR  MATERIAL   TYPE  OF 

707 

C 

LITHOLOGY 

708 

DIMENSION  MAJTS  (10) 

709 

COMMON    IWELL(IOOO) » IVAR ( 8 ) • TV ALU ( 8 ) . IADV ( 8 ) .WVAR ( 1000 ) » ISTOR ( 21 ) 

710 

COMMON  HEAD(2»252).TITLE(2»108).JWELL(1000)»IADVR(8)»STRAT(2Q 

.8  ) 

711 

COMMON   I  SEC  » I  SEC  1 • I SEC2 

712 

100 

FORMAT (  1    '»    'HOLE  NO. '♦ 2X      LAYER   NO '  »2X ♦ 1  MAJOR   TYPE'»3X.    'MINOR  TY 

713 

1  nC  1       CvtClCwn      A  Y  .  1  ftCOTU  1   .  5  V       1  TUT  ^1/  diTC  C  1       cv       IMCT      Turn/  (ma     1  r\n  \ 

1 P  t  1  »  5a  » 1  £  L  t  V  1  » *♦  a  »  'DtPTn1  Oa»  '  TnlCKNcSS'  »5X»  '  N  t  T    THICK ( MAJOR) 

NET 

714 

TTu  trk1  IMlMADI  I  \ 
t  1  n  1  LIn  1  rl  1  'NUK  J  I 

715 

101 

FORMAT!  »2A»I3'BA»I2»lUA»I5»10A»r6tl»2X»r6il»6A»F6«2) 

716 

102 

FORMA)  (  1    1  ifXiPtoXiU  i&cXf  l7tOAtroil92XtFoii  >6Xt  F6i2  ) 

Tit 

717 

103 

718 

104 

719 

Lin  t  T  c  /  a»i  no  i 

720 

C 

T  f"  AM  P\      f  C     TUP     Ml  IMRFD     OP     M  A  TFD  I  Al       T  V  D  P  Q     cCi   FT  TCP* 

I^UINU    1  o     Int.    INUnotK    UP    nA  1  LK1  ML     1  Trto    olLl^.  1  cUi 

721 

UU    O    Is— iilLUINU 

722 

C 

C  P  T  C        iTCcT      Cni  iAI  T '  (~\     UATFD  T  Al        ^~  Pi  P>  P      cFl    CrTCTH      /  n      p\  T  f  T  T  \ 

723 

J  1  CO  1  -  nA  J  1  Oils/ 

724 

c 

cf  arth   tmroi  icm   f  Aru    cpi  FnFn   tf<;t   hDI  f 

725 

726 

c 

72  7 

TAI  1      CFPfLi     (mi  A  Y  ■  f»  Q  F  1  V/  ■  T  T  ) 

72  8 

T  CFT  ?  -  I  I«IF  1   1    t  J  r  \ 
1  OC^t-iV\tLL  I  Is.  J 

72  9 

c 

DpAH     1    A  VFP     riAT  A     CpAM     n  T             PHP     FATW     uHl  P 
KCMU     LA  T  t  r\     UA  1  A     r  iMJi*I    UlOL     rUrv     C.Mv.n  PlvLt 

730 

731 

732 

J=l 

733 

SUM2  =  0 . 0 

734 

SUM3  =  0  *  0 

735 

5  2 

I WLL= I  WELL ( I C ) 

736 

c 

SEARCH  ALL   LAYERS   IN   EACH  HOLE 

737 

53 

DO   U    J  J  =  1  t  NL  AY 

738 

c 

SET    TYPE   ECUAL    TO  MATERIAL   CODE   FOR  EACH   LAYER    (5  DIGIT) 

739 

50 

TYPF=STRAT ( J.3 ) 

cum    tc    THF    TOTAL    TH  TTK^F  sc;    TNrLUuTMft    THF    LAYER    RFTNft  ^FARTHFH 

74  1 

SUM=SUM+STRAT(J»2) 

7b.? 

r 

SETS   KTEST   TO  A   TWO  DIGIT   CODE  OF  MAJOR  MATERIAL  TYPE 

KTEST=IFIX(TYPF) /1000 

"7  A  A 

1- 

THTTK     I  t;    THF    THT  f^MFs^    HP    IAYPR    R  F  T  M  f~i  sFAPfHFn 

745 

Iniv-ls-olKAl  I  J  »  a  J 

746 

c 

1  CO  1     lr    MAJUK    MAItKiAL     lYrt     lo    IN    ocLCL  1  tU    LAY tK 

747 

I  U 1  =  1  T re /  i  UUUU • 

748 

749 

30 

ITYPh=TYPt/1000» 

750 

IDS- i  1 irt     1 i 1 Trt' 1U*1U 1 

751 

Ir II Ut*i 1  ill Jl i Jt 

1  c  0 

752 

31 

T  T  C  C  T  —  T  Pi  0  .*.  1  H 

753 

1000 

I  F ( I  TEST— J TEST  J  32 ♦ 1 »32 

754 

32 

I  1  TP  t -IYPt/10. 

75  5 

lUt-1  1  Trt"l  I  1  I rt / 10*1 U ) 

756 

1001 

IrlID4— ^)35»33»J5 

757 

33 

I  TEST  = I D4+ 10 

758 

TC/TTCCT        ITCf  T  1  /■      1  /< 

I  r  1  ITt5T—JIE5T)H»i»'i 

759 

34 

1  1  YHt=  1  YPt/lOU. 

760 

1UZ3=1  IYPb  —  I  I  T  T  r  t  /  1  (J  0  *  1  0  U  ) 

76 1 

Ir( IDZ3- JTtSTJlOUltlilQOl 

762 

35 

I  T YPE  =  T YPE 

763 

ITEST=ITYPE-(ITYPE/100*100) 

764 

IF( ITEST-JTEST 14*1,4 

765 

1002 

IF(KTEST-JTEST)2»1.2 

766 

c 

SUM1    IS   THE  DEPTH  TO   TOP   OF   LAYER   WHERE  MAJOR   TYPE    IS  FOUND 

767 

1 

SUM1=SUM-STRAT ( J • 2 ) 

768 

c 

SUM2    IS   THE   TOTAL  THICKNESS  OF  LAYERS  WITH  MAJOR  TYPE 

769 

SUM2=SUM2+THICK 

770 

c 

ELEV    IS   ELEVATION  OF  TOP   OF   LAYER   CONTAINING  MAJOR  TYPE 

771 

ELEV=GSELV-SUM1 

772 

c 

SETS   TYPE    TO  AN    INTEGER  NUMBER 

773 

I  TYPE  =  TYPE 

774 

c 

LABELS   HOLE  N0.»    LAYER   N0.»   MAJOR   MATERIAL   TYPE •    ELEV.   UF  TOP 

OF 

775 

C 

LAYER »    DEPTH  TO  LAYER  »    LAYER  THICKNESS 

776 

WRITEI6.101)    IWLL»           JJ  » I  TYPE i ELEV  »SUM1 » THICK 

777 

c 

TEST   FOR  PUNCHED  OUTPUT 

778 

IF( IPUNH14.4.7 

779 

c 

PUNCHES  HOLE  NO.»   LAYER  NO..    MAJOR  MATERIAL   TYPE.   ELEV.   OF   TOP  OF 

780 

c 

LAYER.    DEPTH  TO  LAYER ♦   LAY  THICK 

781 

7 

WRITE(7.8)IWLL»JJ»I TYPE . ELE V . SUM  1 » TH I CK 

782 

8 

FORMAT ( • MAJOR '.3I6.2F8.1.F10.2) 

783 

GO  TO  4 

784 

c 

BREAK  DOWN   TYPE  OF  MATERIAL  CODE   TO   TWO  DIGITS 

785 

2 

INUM=TYPE/1000. 

786 

I TYPE=TYPE 

787 

c 

BREAK  DOWN  TYPE  OF  MATERIAL  CODE   TO  MINOR   TYPE  ONLY 

788 

c 

TEST    IF  MINOR   TYPE  OF  MATERIAL  CODE  CALLED   IS   IN  LAYER 

789 

INUM=TYPE-(FLOAT ( INUM)*1000. ) 

790 

IADJ=INUM/100 

791 

I  F ( I ADJ-1 ) 21 .20 .21 

792 

20 

KTEST=INUM/10 

793 

GO  TO  22 

794 

2  1 

KTEST=INUM-( IADJ*100) 

795 

22 

IF(KTEST-JTEST)4»3»4 

796 

c 

SUM1    IS  THE  DEPTH  TO  TOP  OF  LAYER  WHERE  MINOR  TYPE  CODE    IS  FOUND 

797 

3 

SUM1=SUM-STRAT(J»2) 

798 

SUM3=SUM3+THICK 

799 

c 

ELEV    IS   THE   ELEVATION   OF   TOP  OF  LAYER  WHERE  MINOR  TYPE    IS  FOUND 

800 

ELEV=GSELV-SUM1 

801 

c 

LABELS  WELL  NO..   LAYER  NO..   MINOR   MATERIAL   TYPE.   ELEV  OF   TOP  OF 

802 

c 

LAYER.   LAYER  THICK 

803 

WRITE(6»102)    IWLL.            J J  .  I  T YPE . E LEV . SUM  1 . TH I CK 

804 

I F ( I PUNH ) 4.4.9 

805 

9 

WRITE(7.10)IWLL»JJ. I T YPE. ELEV. SUM  1. THICK 

806 

10 

FORMAT (  'MINOR*  .3I6.2F8.1.F10.2) 

807 

4 

J  =  J+1 

808 

c 

TEST    IF  NET   THICKNESS  OF   LAYERS  WITH  MAJOR  MATERIAL    IS  SELECTED 

809 

IF(SUM2)55»55»54 

810 

c 

WRITE  NET   THICKNESS  OF  LAYERS  WITH  MAJOR  MATERIAL  TYPE 

811 

54 

WR I TE ( 6  » 103 ) SUM2 

812 

55 

IF ( SUM3 15.5.56 

813 

56 

WRITE(6.104)SUM3 

814 

5 

CONT INUE 

815 

6 

CONT I NUE 

816 

RETURN 

817 

END 

818 

c 

819 

c 

820 

SUBROUTINE  TAPES 

821 

c 

SUBROUTINE   TAPES  READS   TAPES  USED    IN  STORAGE  AND  RETRIEVAL  PROGRAM 

822 

c 

TAPE  NUMBER   3   FOR  WELL  VARIABLE  DATA. 

823 

c 

TAPE   DRIVE  NUMBER  9  FOR   LAYER  DATA. 

824 

COMMON    IWELL(IOOO) »IVAR(8)  »TVALU(8)  »IADV(8)  .WVARI1000)  .ISTOR121) 

825 

COMMON  HEAD(2.252) »TITLE(2.108) . JWELL  t 1000)  »IADVR(8)  .STRAT120.8) 

826 

COMMON   I  SEC  .  I  SEC  1 . 1 SEC2 

827 

ISEC1=1 

828 

DO  10  1=1.21 

829 

c 

READS   TEST  HOLE  DATA  VARIABLES  FROM  TAPE 

830 

READ(8.END=11)WVAR 

831 

c 

WRITES   TEST  HOLE  DATA  ON   DISC  FILE 

832 

10 

WR I TE I  2 '  I SEC1 ) WVAR 

833 

1 1 

DC  20  J=1.1000 

834 

c 

READS   LAYER  DATA  VARIABLES   FROM  TAPE 

835 

READ (9.END  =  30) I M  » I . STRAT 

836 

c 

SELECTS  LAYER  VARIABLES  FOR  SPECIFIED  TEST  HOLES 

837 

IF( IM130.30.20 

8*8 

c 

WR I TE ( LAYER  DATA  ON  DISC) 

839 

20 

WRITE( 1 ' I M ) STRAT 

840 

30 

RETURN 

841 

END 

842 

c 

843 

c 

844 

c 

SUBROUTINE   JULDY  COMPUTES  JULIAN   DATE   FROM  CALENDAR  DATE   AND  MILITARY 

845 

c 

TIME. 

40 


646 

SUBROUTINE   JULDY    ( I  MO • I DY  .  I YR » I HR » T I  ME ) 

847 

NYR=IYR/4 

848 

NYR=NYR*4 

849 

IF    (NYR-IYR)  2.1.2 

850 

1 

IB  =  1 

851 

GO  TO  3 

852 

2 

IB"2 

853 

3 

JDAY»0 

854 

DO  10   1=1  » I  MO 

855 

GOTO   I4.5.6.5.9.5. 9.5.5.9. 5.9) »I 

856 

4 

MO»0 

857 

GO  TO  10 

858 

5 

MO  =  31 

859 

GO  TO  10 

860 

6 

GO  TO   ( 7.8 ) » IB 

8  CI 

7 

MO  =  29 

862 

GO  TO  10 

863 

8 

MO  =  28 

864 

GO  TO  10 

865 

9 

MO  =  30 

866 

10 

JDAY=JDAY+MO 

867 

JDAY* JDAY+ I DY 

868 

DAY= JDAY-1 

869 

DAY=DAY*1440. 

870 

HR*(  IHR/100)*60 

871 

HRl=IHR-( ( I HR/ 100)* 100 ) 

872 

TIME=DAY+HR+HR1 

873 

102 

RETURN 

874 

END 

875 

C 

876 

C 

877 

C 

SUBROUTINE  WATLV  COMPUTES  GROUND  WATER   ELEVATION  AT  SELECTED  TEST  HOLES 

878 

C 

WHEN  DATA  FROM  THE  SATURATED  ZONE    IS  RETRIEVED 

879 

SUBROUT INE  WATLV (NWELLS . I  WELL .WAT EL ) 

880 

C 

WATER  LEVEL  PROGRAM  COMPUTES  AVERAGE  GROUND  WATER  ELEVATION   FOR  A  SELECTED 

881 

c 

TIME  PERIOD. 

882 

DIMENSION    I  WELL ( 1000 ) .WATEL (1000 ) 

883 

REWIND  4 

884 

JWELL1=0 

885 

c 

BEGINNING    (IMOt    IDY.    IYR.    ITIM)    AND  ENDING    (JMO.   JDY.    JYR.   JTIM)    MONTH » 

886 

c 

DAY.    YEAR.   MILITARY   TIME   OF  PERIOD  FOR   WHICH  AVERAGE   WATER   LEVEL  IS 

887 

c 

COMPUTED   IS  READ  FROM  DATA  CARD 

888 

1 

READ(5.101 ) I MO. IDY. IYR. ITIM. JMO. JDY. JYR .JTIM 

889 

101 

FORMAT ( 1615) 

890 

WRITE ( 6 . 110 ) I MO. IDY. I YR . I T I M . JMO . JDY . J YR . J T I M 

891 

110 

FORMAT (  ' 0  *  .  25 (  1  * ' ) » *  T IME   PERIOD  FOR  SEASONAL    '  .  I  2  .  1  /  '  .  I  2  .  1  /  •  .  I  2  .  1  / 

892 

1 

'.14.'    TO    '  »  12  .  ' / 1  .  12  .  1  /  •  .  I  2  .  '  /  1  »  14.  IX  .25 (  ' *'  ) .// ) 

893 

1 1 YR  = I YR 

894 

IF( IMO) 107.107.5 

895 

c 

SUBROUTINE   JULDY  COMPUTES   TIME    IN  MINUTES  FROM  MONTH.    DAY.   YEAR.  AND 

896 

c 

MILITARY  TIME 

697 

c 

BT IM-BEG I NN I NG   TIME    (JULIAN   DATE    IN  MINUTES)    OF   PERIOD  FOR  WHICH  WATER 

898 

c 

LEVELS  ARE  AVERAGED 

899 

5 

CALL   JULDY ( IMO.IDY.IYR.ITIM.BTIM) 

900 

IF( IYR-JYR12.3.3 

901 

c 

EYEAR-WHEN  PERIOD  FOR   WHICH  WATER   LEVELS  ARE   AVERAGED  CROSSES  END  OF  YEAR 

902 

c 

THEN  EYEAR    ( END  OF   FIRST  YEAR.   JULIAN  DATE    IN  MINUTES)    IS  CALUCLATED 

903 

2 

CALL  JULDYl 12.31 »IYR. 2359. EYEAR) 

904 

500 

I  I YR  =  I IYR  +  1 

905 

IF( IIYR-JYRJ501.502.502 

906 

c 

STIME-WHEN   PERIOD  FOR   WHICH  WATER   LEVELS  ARE  AVERAGED   SPANS  A  FULL  YEAR 

907 

c 

THAT   YEARS  TIME    (JULIAN   DATE    IN  MINUTES)    IS  CALCULATED 

908 

501 

CALL   JULDY (12. 31.11 YR .2359. STI ME) 

909 

EYEAR=EYEAR+ST IME 

910 

GO  TO  500 

911 

c 

EET I M-THE   JULIAN  DATE    IN   MINUTES  FROM  JANURARY    1   OF   THE   FINAL  YEAR    TO  THE 

912 

c 

END  OF   THE  SELECTED  TIME  PERIOD 

913 

502 

CALL   JULDY ( JMO . JDY . JYR . JT I M . E ET I M ) 

914 

ET IME  =  EET  IM 

915 

ETIM=EETIM+EYEAR 

916 

GO  TO  6 

41 


917  C  ETIM-END  OF  SELECTED  PERIOD    (JULIAN  DATE    IN  MINUTES) 

918  3  CALL   JULDY( JMO  •  JDY . JYR » J T I M . E T I M ) 

919  ETIME=ETIM 

920  C  T T I M-TOTAL   TIME    (MINUTES)    OF  SELECTED   TIKE  PERIOD 

921  6  TTIM=ETIM-BTIM 

922  C  RATE   OF  WATER   LEVEL  CHANGE    IS  COMPUTED  FROM  ELEVATIONS  READ   FROM  TAPE 

923  C  WITHIN   THE  SELECTED  TIME  PERIOD. 

924  DO  1110  I=1.NWELLS 

92  5  IF ( I  WELL (  I  1-JWELL1 17.1000, 8 

926  7   REWIND  4 

927  JWELL1=0 

928  8   KKK  = IWELL (  I  ) 

929  JBRN=1 

930  JBN=1 

931  STIM=0.0 

932  SGH=0.0 

933  TIMZ=0.0 

934  TM=0.0 

935  10   READ(4»END=107)    NWE LL ♦ NCODE . N YR • NMO » ND Y » NT  I M » ELEV 

936  IFINWELL-IWELLU  ))10»20»85 

937  20  CALL  JULDY ( NMO » NDY . NYR . MT I M , T I  ME ) 

938  GO  T0.(  3000.3001  )  .JERN 

939  3000  JBRN=2 

940  IF(  IYR-NYP.)86»30»3010 

941  3001    IF( IYR-NYRI40.30.3010 

942  30  IF(BTIM-TIME)3003,3002»3010 

943  3002  WGH=0.0 

944  ELEV2=ELEV 

945  TIMEB=TIME 

946  GO  TO  3011 

947  3003  IF(NYR1-NYR)3005»3004»3004 

948  3004  JBN=2 

949  GO  TO  40 

950  3005  JBN=3 

951  GO  TO  40 

952  3010  NYR1 =NYR 

953  GO  TO  10 

954  40   BACKSPACE  4 

955  BACKSPACE  4 

95  6  READ (4 »END= 107 )N WELL t NCODE* NY R.NMOtNDY  »NTIM» ELEV 1 

957  CALL  JULDY(NMO.NDY.NYR.NTIM.TIMA) 

958  45  READ(4.END=107)KWELL.KCODE.KYR.KMO.KDY.KTIM.ELEV2 

959  IFUWELL-IWELL  (  I  )  )86»46»87 

960  46    IF( (KYR-NYR)-l )3014. 3013. 3012 

961  C  A  BREAK    IN  RECORD  GREATER   THAN  ONE  YEAR   IS  PRINTED   IN  OUTPUT 

962  3012  WRITE(6»109)IWELL(I)  »NMO » NDY » NYR » NT  I M » KMO . KDY  ,  K YR  ♦  KT I M 

963  109   FORMAT (  '    WELL  NUMBER    •  .14,'    BREAK    IN  RECORD  FROM    • »3I3iI5»'  TO 

964  1313.15) 

965  GO  TO  1005 

966  3013    IF(KMO+( 1 2-NMO ) - 12 ) 30 14  .  3012  .3012 

967  3014  CALL   JULD Y ( KMO  .  KDY , KYR . K T I M . T I MEB ) 

968  GO  T01206. 207. 3006) .J3N 

969  3006  CALL   JULDY ( 1 2 » 3 1  .NYR » 2 3 59  . TM ) 

9  70  D  I  FF= (  ( TM-TIMA )+BT IM )  +  ( T I MEB-BT I M ) 

971  DIFF1=TIMEB-BTIM 

972  GO  TO  206 

973  206  CALL   JULDY( 12.31  .NYR. 2359. TM) 

974  DIFF=(TM-TIMA)+TIMEB 

975  DIFF1=(TM-BTIM)+TIMEB 

976  GO  TO  208 

977  207  DI FF  =  T I MEB-T I  MA 

978  DIFF1=T IMEB-BT IM 

979  208   SLOPE= ! ELCV2-ELEV1 ) /DI FF 

980  GHA=ELCV2-(SL0PE*DIFF1 ) 

981  D I FF=D I FF 1 

982  STIM=STIM+DIFF 

983  IF(STIM-TTIM) 1002. 1011. 1010 

984  1010  DIFF=DIFF-(TM-ETIM+TIMEB) 

9^5  GHB=ELEV2-SL0PE* (TM-ET IM  +  T I MEB  > 

986  GO  TO  1012 

987  1011  -GHB=ELEV2-SL0PE*  (  T  IMEB-ET-IM! 
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1000 

DIFF=( (TM-TIMZ1+TIMX) 

1001 

GO  TO  203 

1002 

201 

DIFF=T IMX-TIMZ 

1003 

203 

STIM=STIM+DIFF 

1004 
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1018 
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1019 
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1  0  5  1 

TF  /  MV'Bl  »mvd  \  one  ,  one;  »orii 

1  V  £  J 

a  n  l 

O  U  A 

1024 

DIFF= ( TMT -TIMCA)+TTMFR 

1025 

^L0PF= ( ELFV2-FLEV1 ) /DI  FF 

1026 

IF(JYR-NYR ) 803  »803»802 

1027 

802 

GHB= FIEV2- ( SLOPE* ( T I MEB-E T I  ME )  ) 

1028 

WGH=( ( ( (GHB+GHA)*»5)*( (TM1— TIMEA)    +FTIMF) l/TTIM) 

1029 

GO   TO  9001 

1030 

803 

GHB*ELEV1+ ( SLOPE* ( ET IME-T I ME A ) ) 

1031 

GO  TO  9000 

1032 

805 

SL0PE=1ELEV2-ELEV1)/(TIMEB-TIMEA) 

1033 

RHR:FLEV>ISI  OPE*ITlMFR-FTIWFI  1  ■ 

•J    'lJ  I —  L  L.  V  C.        *   J  L  V  r   I,  "    \    |    i  i"l  L  U       L.    1    1  I'lL   /  1 

10  34 

TF ( FTTMF-T  IMFA  1Q000 .QOOQ «Q000 

1035 

9009 

FT  IME  =  T IMEB 

L.    I    i  rIL        1    1  I  I  L  U 

1036 

9000 

WGH=(  (  (  ( GHB+GHA ) *• 5 ) * ( ET I  ME— T I ME A )  J/TTIM) 

1037 

9001 

WGHP=SGH+WGH 

10  3  8 

GO   TO    1 006 

1  0  7  9 

A  ^  J  7 

r 

(n/MFM     MO    PFTnOPl              FHllMH     CAP     CCI  prrcn     T  IMP     DPD  IOH     TWTC     TC     DPTMTFH     TM  OIITDllT 

1 040 

1041 

103 

FORMAT!'           WELL  N0«    '»I4»'    NO  RECORD  FOR  PERIOD  OF    ' »3I3»I5t'  TO 

1042 

1    '  »  3  I  3  »  I  5  ) 

1043 

GO  TO  100  5 

1 044 

WHFN    <iFLFC"TFr)    T I MF    PERIOD    l>TAHTe;    RFFORF    RFrORD^    ON    TAPE    THTC,    I  c;    PRINTED  IN 

1  0^*5 

A  w  "t  »1 

c 

OUTPUT. 

1UHO 

86 

WRITE(6tl04) IWELL(  I  ) 

1 U  4  ' 

104 

FORMAT ( '    't'   WELL  MO.    'tI4.'    PERIOD  OF  RECORD  STARTS  BEFORE  ACTUAL 

1048 

1REC0RD ' ) 

GO  TO  1005 

1  (Tin 

C 

WHEN  SELECTED  TIME  PERIOD  GOES  BEYOND  RECORDS  ON  TAPE   THIS   IS  PRINTED  IN 

1  0  S  1 

c 

OUTPUT 

1  V7  £ 

87 

WRITE(6.105) IWELL( I ) 

1  n  e.  *a 

105 

FORMAT ( 1    't'    WELL  NO.    '»I4»'    PERIOD  OF   RECORD  BEYOND  ACTUAL  RECORD 

i  o  s^a 

1  '  ) 

10  5  5 

C 

NUMBER   99999  STORED    IN  WATEL    UKK)    INDICATES  A  MISSING  RECORD  FOR  THESE 

1056 

C 

SELECTED  TIME  PERIODS 

1057 

1005 

WGHP=99999. 

1058 

1006 

WATEL(KKK)=WGHP 
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1059  1000   JWELLl= IWELL ( I ) 

1060  1110  CONTINUE 

1061  107  RETURN 

1062  END 

1063  C 

1064  C 

1065  C  END  OF  RETRIEVAL  PROGRAMS,   STORAGE  PROGRAMS  FOLLOW. 

1066  C 

1067  C  ****  STORAGE   PROGRAM   FOR   TEST  HOLE  DATA 
1066                   DEFINE   FILE   1(21  1 1000  »U  » I M  ) 

1069  DIMENSION   WELLl < 1000 ) »WELL ( 4 ) .NN ( 4 ) 

1070  DIMENSION    I  START (21) 

1071  DO  105  1=1.21 

1072  105  ISTART(I)=0 

1073  C  READ  FROM  CARD  THE  NUMBER  OF   TEST  HOLE  VARIABLES  TO  BE  STORED  ON  DISC 

1074  READI5.100)  NOVAR 

1075  100   FORMAT ( 14 ) 

1076  C  READ  A  VARIABLE  NUMBER*    FOUR   VARIABLE  VALUES.   AND  FOUR  TEST  HOLE 

1077  C  NUMBERS  FROM  CARDS 

1078  99  READ( 5  .200  INVAR  .  (NN < I  )  » I  =  1  .4  )  .  < WELL (  I >  » 1  =  1 .4) 

1079  200   FORMAT!  I  2 » 4  I  5 • 4F 1 1 . 0 ) 

1080  IF(NVAR)4.4,1 
1031                1    I M=N VAR 

1082  IF( ISTART(NVAR).EQ.O)    GO   TO  6 

1083  READ ( 1 ' IM )  WELLl 

1084  IM=IM-1 

1085  GO  TO  8 

1086  6  DO  7  1=1.1000 

1087  7  WELLl ( I ) =0. 

1088  8  DO  3  1=1.4 

1089  IF(NN(  I  )  13.3.2 

1090  2  IW=NN(I) 

1091  WELL  1 ( I W ) =WELL (  I  ) 

1092  3  CONTINUE 

1093  C  WRITE  TEST   HOLE  VARIABLE  ON  DISC 

1094  WRITE(l'IM)  WELLl 

1095  ISTART(NVAR)=1 

1096  GO  TO  99 

1097  C  WRITE  TEST  HOLE  VARIABLE   DATA  FROM  DISC   TO  TAPE 

1098  4  IM=1 

1099  DO  5    1=1. NOVAR 

1100  READ(l'IM)  WELLl 

1101  WRITE(9)  WELLl 

1102  C  WRITE  VARIABLE  NUMBER .    TEST  HOLE   NUMBER.   AND  VARIABLE  VALUE'ON  PRINTER 

1103  WR ITE (6. 101)1. ( J  .WELLl ( J ) . J  =  l . 100  0 ) 

1104  101   FORMAT (' 0 VARIABLE  NUMBER '  » I  2 » / / » 7 ( I  4  .  2X  .  F 1 1 . 0 ) ) 

1105  5  CONTINUE 

1106  WRITE16.102) 

1107  102   FORMAT    I'l'.'TEST  HOLE  VARIABLE  DATA  HAS  BEEN  STORED  ON  TAPE') 

1108  STOP 

1109  END 

1110  C 

1111  C  ; 

1112  C 

1113  C  ##♦*   STORAGE   PROGRAM  FOR   LAYER  DATA  **** 

1114  DEFINE   FILE    1  ( 1 000 , 16 1  . U . I M ) 

1115  DIMENSION   STRAT ( 20 .8 )  .NWELLS I  1000  ) 

1116  C  NOWEL    IS   THE  NUMBER  OF   TEST   HOLES  FOR  WHICH  DATA    IS  STORED 

1117  READ ( 5 . IOC ) NOWEL 

1118  100   FORMAT ( I  4 ) 

1119  ICOUN=0 

1120  DO  3    IC=1. NOWEL 

1121  C  NWELL    IS   THE   TEST  HOLE  NUMBER 

1122  C  STRAT    (l.J)    IS  THE  LAYER  VARIABLES  STORED   IN  EACH  LAYER 

1123  C  NOLAY    IS   THE  NUMBER  OF   THE   LAYER   FOR  WHICH  VARIABLES  ARE  STORED. 

1124  READ (5.101) NWELL  .(STRAT(l.J)  »J=1.8) .NOLAY 

1125  101   FORMAT! I4.2F4.0.1X.F5.0.2F3.0.2F7.0.F10.6.I4) 

1126  IF(NWELL)200.200.300 

1127  300  ICOUN=ICOUN+l 

1128  IFINOLAY-U3.2.1 

1129  1   READ! 5 » 106 )  (  (STRATI  I . J)  . J  =  l .8 )  .  1  =  2. NOLAY) 


44 


1130  106  FORMAT(4X.2F4.0»lX,F5.0.2F3.O»2F7.0»Fl0.6) 

1131  2  WRITE! 1 'NWELL ) NOLAY .STRAT 

1132  NWELLSl  IO=NWELL 

1133  3  CONTINUE 

1134  C  SEQUENCE   THE  WELLS   IN  ASCENDING  ORDER 

1135  200  NUM=ICOUN-l 

1136  DO  5   J=l » ICOUN 

1137  DO  5  I=1»NUM 

1138  IF(NWELLS( I )-NWELLS( 1+1 ) ) 5. 5  .4- 

1139  4  HOLD=NWELLS< I ) 

1140  NWELLSl I ) =NWELLS ( 1+1 > 

1141  NWELLSl 1  +  1 )=HOLD 

1142  5  CONTINUE 

1143  C  WRITF  LAYER  DATA  ON  TAPE  AND  PRINTER 

1144  DO  7  IC=l»ICOUN 

1145  I M  =  NWELLS I  I C  > 

1146  READ ( 1 ' IM ) NOLAY »STRAT 

1147  IU=IM-1 

1148  WRITE(8)    I U  »  NOLA Y  t  STR AT 

1149  DO  6  I=l»NOLAY 

1151  103  FORMAT ( •    '  .  2  I  5 » 8 F10 . 1 ) 

1152  7  CONTINUE 

1153  IM=0 

1154  WRITEI8) IM.NOLAY.STRAT 

1155  WRITE(6»104) 

1156  104  FORMAT (• 1 '»' LAYER  DATA  STORED  ON   TAPE  1  ) 

1157  STOP 

1158  END 

1159  C 

1160  C 

1161  C  ****  WATER  LEVEL  DATA  STORAGE  PROGRAM  **** 

1162  C  CALLS  SUBROUTINE  (TESTED 

1163  C  PROGRAM  STORES  PUNCH  CARD  TAPE  DOWN  DATA  AS  WATER  LEVEL  ELEVATIONS  ON 

1164  C  MAGNETIC   TAPE • 

1165  C  MASTER  TAPE  OF  WATER  LEVEL  ELEVATIONS   IS  CONSTRUCTED  USING  PIPE  ELEVATIONS 

1166  C  AND  TAPE  DOWN  DATA  STORED  ON  CARDS  OR  TAPES. 

1167  C  WATER  LEVEL  DATA  FOR  SELECTED  TEST  HOLES  CAN  BE  UPDATED  USING  THIS 

1168  C  PROGRAM • 

1169  C  TAPE  PARAMETERS  READ   IN  BY  DATA  STATEMENT 

1170  C  WORDS-NUMBER  OF  WORDS  WRITTEN  PER  RECORD 

1171  C  LRECL-LENGTH  OF  RECORD  SIZE 

1172  C  BLKSIZ-BLOCK  SIZE  PARAMETER 

1173  C  BPI-BITS  PER   INCH  FOR  OUTPUT  TAPE 

1174  C  RGAP-RECORD  GAPE  BETWEEN  RECORDS  (INCHES) 

1175  REAL*4  LRECL ♦ I NPBK 

1176  DIMENSION  KWL ( 4 ) »ELVX ( 4 ) t EL VN ( 100 0 ) » I T EST ( 1 000 ) » JWELL ( 3 ) 

1177  DIMENSION    I WN ( 100 ) , IMO ( 1 0 0 ) ♦ I DY ( 1 00 ) . I YR ( 1 00 ) » ELV ( 100 ) 

1178  DIMENSION    I  TAPE ( 3 ) »NMO ( 3 ) . NDY ( 3 ) » NYR ( 3 ) » NT  I M ( 3 ) • EL  I  3 ) 

1179  COMMON  NYR»NMO»NDY»NTIM»EL» IYR» I  MO ♦ I DY » ELV , ELVN , I TES T ♦ KT APE . 

1180  *NCODE»ICOUNT.IPAGE.IWN 

1181  DATA  WORDS»LRECL*BLKSIZ»BPI »RGAP/8. »24. »2400. »1600. t .75/ 

1182  DATA  JWELL/3*0/ 

1183  DO  1  1=1,1000 

1184  I  TEST ( I )  =  1 

1185  1  ELVN ( I  )  =0.0 

1186  C  READ   IN   TOP  OF  PIPE  ELEVATIONS  FOR  ALL   TEST  HOLES. 

1187  C  IV-LAST  CARD  TEST 

1188  C  KWL   (I)-TEST  HOLE  NUMBERS. 

1189  C  ELVX    (I)-TOP  OF  PIPE  ELEVATIONS  AT  EACH   TEST  HOLE. 

1190  5  READ ( 5  » 1000 ) I V  » ( KWL ( I >  1 1  =  1 »4 ) »  t ELVX ( I ) » I  =  1 • 4 ) 

1191  1000  FORMAT ( 12 ,415 »4F11.2 ) 

1192  IF(IV.LE.O)   GO  TO  15 

1193  DO  10  1=1.4 

1194  IF(KWL( I ) .LE.O)   GO  TO  10 

1195  KKK=KWL ( I ) 

1196  ELVN(KKK)=ELVX( I ) 

1197  10  CONTINUE 

1198  GO  TO  5 

1199  C  READ   IN  PIPE  ELEVATIONS  CHANGES 

1200  15  1=1 

1201  C  IWN   (I)-TEST  HOLE  NUMBER 
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1202  C  I  MO    (I).    IDY    (I)t    IYR    (I)    IS  THE  MONTH »    DAY  »   AND  YEAR  OF   PIPE  ELEVATION 

1203  C  (ELV    (I))  CHANGES. 

1204  16  READ ( 5  » 1 00 1 ) I WN ( I ) » I MO ( I ) • I DY ( I ) » I YR ( I ) » ELV ( I ) 

1205  1001  F0RMAT(4I5»F10.2) 

1206  IF ( IWN (  I  )  .LE.O )    GO   TO  17 

1207  I  TEST ( IWN ( I ) ) =0 

1208  1=1+1 

1209  GO  TO  16 

1210  17  N0BS1=I 

1211  C  NTAPE-NUMBER  OF   TAPE   DRIVE  UNITS  USED. 

1212  C  ITAPE    (I)-TAPE  DRIVE  UNIT  NUMBERS . 

1213  C  <TAPE-CREATED  OR  UPDATED   TAPE  UNIT  NUMBER. 

1214  READ (5»1002)NTAPEt(ITAPE(I),I=l»NTAPE)»KTAPE 

1215  1002  FORMATI5I5) 

1216  C  READ  WELL  NUMBER  FROM   TAPES  FOR    I N I T I L I ZAT I  ON • 

1217  1005  DO  1006  IT=1,NTAPE 

1218  ITP=ITAPE( IT) 

1219  READ( ITP»103) JWELLt IT) 

1220  1006  CONTINUE 

1221  1007  WR I TE ( 6  »9999 ) ( JWELL ( I ) » I  =  1 1 NT APE ) 

1222  -9999   FORMAT (  1 0*****    '  »  3  I  10  » '  *****•/) 

1223  C  FOLLOWING   TEST   STATEMENTS  SELECT   DATA  FROM   TWO  OR   THREE  DATA  SETS  AND 

1224  C  ARRANGES   THEM    IN  ORDER  OF  ASCENDING  TEST  HOLE  NUMBERS   FOR  CREATING  OR 

1225  C  UPDATING  DATA  TAPE. 

1226  IF  I JWELL ( 1 ) .EQ.O.AND. JWELL ( 2 > .EQ.O.AND. JWELL! 3) .EQ.O ) GO   TO  2000 

1227  IF (JWELL ( 1 ) .EQ.O.AND. JWFLL ( 2 ) .EQ.O.AND. JWELLI 3 ) .GT.O ) GO   TO  1060 

1228  IF (JWELLI 1 ) .EQ.O.AND. JWFLL ( 2 ) .GT.O. AND. JWELLI 3) .EQ.O ) GO   TO  1050 

1229  IF  (JWELL  (  1  )  .C-T.O.AND.  JWELL  (  2  )  .EQ.  0. AND.  JWELL  (  3)  .EQ.O  )G0   TO  1030 

1230  IF ( JWELL ( 1 ) .EQ.O.AND. JWELL ( 2 ) .EQ. JWELL ( 3 ) )    GO  TO  1040 

1231  IF(JWFLL(  1  )  .EQ.O.AND. JWELLI 2 ) .LT. JWELLI 3)  )    GO  TO  1050 

1232  IF ( JWELL ( 1 ) .EQ.O.AND. JWELL I  2 >  .GT. JWELL I  3 )  )    GO  TO  1060 

1233  IF(JWELL(2) .EQ.O. AND. JWELL(l) .EQ.JWELLI3) )    GO  TO  1020 

1234  I F I JWELL I  2 ) . EQ « 0 . AND. JWELL I  1 ) .LT. JWELL ( 3 ) )   GO  TO  1030 

1235  IF (JWELL ( 2 ) .EQ.O.AND. JWELL ( 1 ) .GT. JWELL I  3 )  )    GO  TO  1060 

1236  IF  !  JWELL I  3  )  .EQ.O.AND. JWELL I  1 ) .EQ. JWELL I  2 )  )    GO  TO  1010 

1237  IF (JWELL ( 3 ) .EQ.O.AND. JWELL I  1 ) .LT. JWELL I  2 )  )    GO  TO  1030 

1238  IF (JWELL ( 3 > .EQ.O.AND. JWELL ( 1 ) .GT. JWELL ( 2 ) )    GO  TO  1050 

1239  IF ( JWELL ( 1 )  .EO. JWELL( 2 ) .AND. JWELL ( 2 ) .EQ. JWELL ( 3 )  )    GO   TO  1009 

1240  IF  I JWELL (  1  ) .EQ. JWELL( 2 ) .AMD. JWELL ( 2 ) .LT. JWFLLI3 )  )    GO   TO  1010 

1241  IF( JWELL ( 1 ) .LT.JWELL ( 2 ) .AND. JWELL ( 1 ) sEQ. JWELL (3 ) )    GO   TO  1020 

1242  IF  (  JWELL  (  1  )  .LT.  JWELLt  2  )  .AND.  JWELL  (  2  )  .EQ.JWELH3  )  )    GO   TO  1030 

1243  IF ( JWELL ( 1 ) .GT. JWELL( 2 ) .AND. JWELL ( 2 ) .EQ. JWELL I  3 )  )    GO   TO  1040 

1244  IF ( JWELL ( 1 )  .GT. JWELL I  2 ) .AND. JWELL I  2 ) .LT. JWELLI 3 )  )    GO   TO  1050 

1245  IF! JWELLI 1 ) .LT. JWELL I  2 ) .AND. JWELL I  2 ) .LT. JWELL (3  )  )    GO   TO  1030 

1246  IF ( JWELL( 1 1 .LT. JWELLI 2 ) .AND. JWELL (2 ) .GT. JWELL ( 3 ) )    GO   TO  1030 

1247  IF  I JWELLI 1 ) .GT. JWELLI 2 >  .AND. JWELL ( 2 ) .LT. JWELL ( 3 )  )    GO   TO  1050 

1248  IF  I  JWELLI 1  )  .GT. JWELL I  2 )  .AND. JWELL (2 ) .GT. JWELL ( 3 )  )    GO   TO  1060 

1249  WRITE! 6.1008 ) 

1250  1008   FORMAT! '0' » 'THERE  ARE   NO   TEST   COMBINATIONS  FOR   THESE   TEST  HOLES.'/) 

1251  GO  TO  2000 

1252  1009  DO  34  IT=1.NTAPE 

1253  I TP= I  TAPE  I  I  T ) 

1254  IFI ITAPE! IT) .EQ.O)   GO   TO  34 

1255  BACKSPACE  ITP 

1256  NWELL1  =  JWELL( ITP  ) 

1257  20   IF(ITP.EQ.l)   GO  TO  21 

125  8  READ  I  I  TP . 103 » END  =  30 ) NWELL • NCODE » ( NMO I  I )  »  NDY I  I  ) » NYR I  I  )  .NT IM (  I  ) . EL ( I 

1259  1).I=1.3) 

1260  I READ=3 

1261  103  FORMAT ( 14. 12.3 ( 1X.3 12 ♦ I  5  .F7.2  )  ) 

1262  GO  TO  22 

12  63  21   READ (ITP. 104 » END  =  30 ) NWELL . NCODE » NMO ( 1 )  .NDY I  1  )  .NYR  I  1  )  » NT  I M I  1  )  .  EH  1 ) 

1264  104  FORMAT  I  14 .412 » I4.F8.2 ) 

1265  I READ= 1 

1266  22    IFINWFLLl.NE. NWELL)    GO   TO  25 

1267  CALL   TESTEHNWELL.NOBSl  .IREAD) 

1268  NWELLl=NWELL 

1269  GO  TO  20 

1270  25   JWELLt I  TP ) = NWELL 

1271  GO  TO  34 

1272  30  JWELL(ITP)=0 
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1273  I  TAPE ( I T ) =0 

1274  34  CONTINUE 

1275  GO  TO  1007 

1276  1010  DO  45  IT=1.2 

1277  I  TP  =  I  TAPE  (IT) 

1278  I F ( I  TAPE ( I T ) • EQ • 0 )   GO   TO  45 

1279  BACKSPACE  ITP 

1280  MWELLl=JWELL( ITP) 

1281  35    IF(ITP.EQ.l)   GO  TO  36 

1282  READ(ITPtl03»END=41 ) NWELL  » NCODE ♦ ( NMO ( I  )  »  NDY ( I ) ♦ NYR I  I  > • NT  I M ( I  >  » EL  (  I 

1283  1)»I=1»3) 

1284  I READ=3 

1285  GO  TO  37 

1286  36  READ ( ITP i 104 »END=4l ) NwELL , NCODE » NMO ( 1 ) f  NDY ( 1 ) »NYR( l)»NTIM(l)iEL(l) 

1287  I READ= 1 

1288  37    IF ( NWELL 1 »NE  «NWELL )   GO   TO  40 

1289  CALL   TESTEL ( NWELL tNOBSl » I  READ) 

1290  NWELLl=NWELL 

1291  GO  TO  35 

1292  40   JWELL (  I  TP ) =NWELL 

1293  GO  TO  45 

1294  41   JWFLLI ITP)=0 

1295  ITAPE(IT)=0 

1296  45  CONTINUE 

1297  GO  TO  1007 

1298  1020  DO  65   I T=  1 » 2 

1299  GO  TO   <  50»51 )  .  IT 

1300  50  ITP=ITAPE<1) 

1301  IF( ITAPE(l) .EQ.OJGO  TO  65 

1302  GO  TO  52 

1303  51  ITP=ITAPE(3) 

1304  IF( ITAPE(3).EQ.0)GO  TO  65 

1305  52  BACKSPACE  ITP 

1306  NWELLl3 JWELL ( I  TP ) 

1307  55    IF(ITP.EQ.l)   GO  TO  56 

1308  READ  I ITP»103»END  =  61 > NWELL • NCODE ♦ <NMO< I )  » NDY ( I ) t  NYR ( I ) i NT  I M ( I ) »  EL ( I 

1309  1).I=1.3) 

1310  I  READ  =  3 

1311  GO  TO  57 

1312  56  READ ( I  TP » 1 04 » END =6 1 ) NWELL ♦ NCODE t NMO ( 1 )  » N D Y ( 1 )  » N Y R ( 1 )  iNTIMl 1 ) »EL( 1 ) 

1313  I  READ= 1 

1314  57   I F ( NWELL1 «NE« NWELL )   GO   TO  60 

1315  CALL   TESTEL (NWELL* NOBS 1 ♦ I  READ) 

1316  NWELL1=NWELL 

1317  GO  TO  55 

131S  60   JWELL1 ITP)=NWELL 

1319  GO  TO  65 

1320  61  JWELL(ITP)=0 

1321  ITAPE(ITP)=0 

1322  65  CONTINUE 

1323  GO  TO  1007 

1324  1030  ITP=ITAPE<1) 

1325  IF ( ITAPEt 1 ) .EQ.O )   GO  TO  1007 

1326  BACKSPACE    I  TP 

1327  NWELL1  =  JWELL ( I  TP  ) 

1328  70   IF(ITP.EO.l)   GO  TO  71 

13  29  READ!  I  TP . 103 . END  =  80 ) NWELL • NCODE » ( NMO < I  ) »  NDY ( I  ) • NYR  (  I  )  » NT  I M I  I ) f EL ( I 

1330  1  )  .  1  =  1  *3) 

1331  I  READ  =  3 

1332  GO  TO  72 

1333  71   READ  (  I  TP  t 104  .  END  =  eO  )  NWELL  .  NCODE  •  NMO  (  1 )  »NDY(  1)  »NYR  (  l)»NTIM(l)tELU) 

1334  I READ= 1 

1335  72    IFdNWELLl.NE. NWELL)   GO   TO  75 

1336  CALL  TESTEL(NWELL»N0BS1»IREAD) 

1337  NWELL1=NWELL 

1338  GO  TO  70 

1339  75   JWFLLI ITP Is NWELL 

1340  GO  TO  1007 

1341  80  JWELL ( ITP)=0 

1342  ITAPE(1)=0 

1343  GO  TO  1007 
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1344 

1040 

DO  100  IT=li2 

1345 

GO  TO< 85.86) » I T 

1346 

85 

ITP=ITAPE(2) 

1347 

IF( ITAPEI2) .EQ.01GO  TO  100 

1348 

GO  TO  87 

1349 

86 

ITP=ITAPE(3) 

1350 

IF( ITAPEI 3) .EQ.0)GO  TO  100 

1351 

87 

BACKSPACE  ITP 

1352 

NWELL1" JWELL ( I  TP ) 

1353 

90 

IF( ITP.EO. 1 )   GO  TO  91 

1354 

READ ( ITP.103.END=96)NWELL»NCODE» ( NMO ( I ) t 

NDY (  I  ) 

»NYR( I  ) »  NT  I  M (  I  )  » 

EL  (  I 

1355 

1)  »  I  =  1 »  3  ) 

1356 

I  READ  =  3 

1357 

GO  TO  92 

1358 

91 

READ ( ITP» 1 04 ♦ END  =  96  )  NWELL  »NCODE»NMO< 1 ) » NDY ( 1 ) » 

NYRl 1 ) tNTIMl 1 ) tELl 1 ) 

1359 

I READ= 1 

1360 

92 

IF(NWELLl.NE.NWELL)   GO  TO  95 

1361 

CALL   TESTELtNWELLiNOBSl » IREAD) 

1362 

NWELLl3NWELL 

1363 

GO  TO  90 

1364 

95 

JWELLI ITP)=NWELL 

1365 

GO  TO  100 

1366 

96 

JWELLI ITP)«0 

1367 

ITAPE ( I  TP ) *0 

1368 

100 

CONTINUE 

1369 

GO  TO  1007 

1370 

1050 

ITP=ITAPE(2) 

1371 

IF( ITAPE(2)  .EQ.01GO  TO  115 

1372 

BACKSPACE  ITP 

1373 

NWELL 1  =  JWELL { I  TP ) 

1374 

105 

IF(ITP.EQ.l)   GO  TO  106 

1375 

READ ( ITP»103»FND«111)NWELL»NCODE» (NMO( I ) 

»NDY( I 

)  » NYR (  I  )  .NTIMI  I  )  »ELl 

1376 

11  )  »I=1»3) 

1377 

I  READ  =  3 

1378 

GO  TO  107 

1379 

106 

READ  (  ITP.104.END  =  111)NWELL»NCODE»NMOU)  » 

NDY (  1  ) 

»  NYR ( 1  )  »NTIM(  1  )  . 

ELI  1 

1380 

») 

1381 

IREAD=1 

1382 

107 

I F ( N WELL 1 »NE • NWELL )   GO  TO  110 

1383 

CALL   TESTELfNWELL.NOBSl » IREAD) 

1384 

NWELL1=NWELL 

1385 

GO  TO  105 

1386 

110 

JWFLLI ITP)=NWELL 

1387 

GO  TO  115 

1388 

111 

JWELL (  ITP  )*0 

1389 

I  TAPE ( ITP)=0 

1390 

115 

CONTINUE 

1391 

GO  TO  1007 

1392 

1060 

ITP=ITAPE ( 3  ) 

1393 

IF(  ITAPE(  3)  .EQ.O  )GO  TO  130- 

1394 

BACKSPACE  ITP 

1395 

NWELL 1 s JWELL ( I  TP ) 

1396 

120 

IF(ITP.EQ.l)   GO  TO  121 

1397 

READ ( ITP»103»END  =  126)NWFLL»NCCDE»  t NMO (  I ) 

tNDY I  I )  » N Y  R I  I  )  » N  T I M (  I  ) 

•  EL  I 

1398 

11  )  •  I  =»  1  •  3  ) 

1399 

I READ=3 

1400 

GO  TO  122 

1401 

121 

READ (  ITP»104»END  =  126)NWELL.NC0DE»NM0l  1)  » 

NDY ( 1  ) 

»  NYR ( 1 )  *  NT  I M ( 1  )  » 

EL  I  1 

1402 

*> 

1403 

I  READ= 1 

1404 

122 

IF ( NWELL1 «NE» NWELL )   GO  TO  125 

1405 

CALL   TESTEH NWELL tNOBSl » IREAD) 

1406 

NWELL 1 =NWELL 

1407 

GO  TO  120 

1408 

125 

JWELLI I  TP ) =NWELL 

1409 

GO  TO  130 

1410 

126 

JWELL ( ITP)=0 

1411 

I  TAPE ( ITP)=0 

1412 

130 

CONTINUE 

1413 

GO  TO  1007 

1414 

C 

COMPUTES   THE  NUMBER  OF   FEET  OF   TAPE  USED 

■  TC  STORE   WATER  LEVEL 

ELEVATIONS 
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1415  C  AND  THE  NUMBER  OF  RECORDS  STORED. 

1416  2000  RECD=ICOUNT 

1417  RECDS=BLKSIZ/LRECL 

1418  INPBK=<  (WORDS*BLKSIZ ) /BP  I  J+RGAP 

1419  FT=( (RECD/RECDS>*lNPBK) /12. 

1420  WRITE(6t2001)    I  COUNT iKT APE  *  FT 

1421  2001   FORMAT! '0' • '$$$$$      '»I7»'    NUMBER   OF  OBSERVATIONS  WRITTEN  ON   TAPE  '»I2» 

1422  *'      $$$$$' /»6X. 'FEET  USED   TO  STORE   DATA  ON   THIS  TAPE    IS  '»F6.D 

1423  REWIND  KTAPE 

1424  STOP 

1425  END 

1426  C 

1427  C 

1428  C  SUBROUTINE   TESTEL  CALCULATES  WATER   LEVEL   ELEVATIONS  AND  WRITES  ON  UPDATED 

1429  C  TAPE. 

1430  SUBROUTINE   TESTE L ( NWEL L , NOBS1 . I  RE AD ) 

1431  COMMON  NYR  <  3 ) »NMO( 3 ) »NDY ( 3  > ♦ NT IM ( 3 >  » EL <  3 ) ♦ I YR ( 100 ) 1 1  MO ( 1 00 > ♦ I DY < 10 

1432  10)  .ELVUOO)  »ELVN  (  1000)  »  I  TEST  (  1000  )  .  K.TAP  E  .  NCODE  .  I  COU.N  T  .  I  P  AGE  .  IWN(  10 

1433  10) 

1434  DO  28    1  =  1. I  READ 

1435  IF( IREAD.EQ.l )    GO  TO  1000 

1436  IF(NYR( I ) .EQ.O)    GO   TO  28 

1437  IF ( I  TEST ( NWELL ) »EQ.  1 )   GO   TO  27 

1438  DO  25  J'l.NOBSl 

1439  C  TEST   FOR  PIPE   ELEVATION  CHANGE 

1440  I F ( NYR ( I ) • GE • I YR ( J ) .AND.NMO ( I ) »GE • IMO ( J ) . AND.NDY ( I ) .GE« I DY ( J ) .AND* 

1441  1NWELL.EO. IWN ( J ) )    GO  TO  26 

1442  25  CONTINUE 

1443  GO  TO  27 

1444  C  NEW  PIPE  ELEVATION 

1445  26   ELVN(NWELD=ELV(  J) 

1446  ITEST(NWELL)=1 

1447  C  WATER   LEVEL  ELEVATION  CALCULATED. 

1448  27   ELI  I  ) =ELVN < NWELL ) -EL (  I  ) 

1449  C  OUTPUT   WRITTEN  ON  UPDATED  TAPE 

145  0         1000  WR I TE( KTAPE* 104) NWELL »NCODE«NMO( I ) » NDY (  I ) t NYR ( I ) »NT I M ( I ) »EL ( I  I 

1451  104  FORMAT( I4.4I2» I4.F8.2 ) 

1452  28  CONTINUE 

1453  RETURN 

1454  END 

1455  C  END  APPENDIX  C 
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APPENDIX  D.-SOURCE  LISTING  FOR  HYDRAULIC-COEFFICIENT 

RANGE-ASSIGNMENT  PROGRAM 


1 

1 
2 

L 
C 

MrrtiNUlA  U 

3 
4 

c 
c 

HYDRAULIC  COEFFICIENT  RANGE  ASSIGNMENT  PROGRAM 

5 

c 

***THIS  PROGRAM  ASSIGNS  HYDRAULIC  COEFFICIENT  RANGE   TO  EACH  LAYER*** 

6 

c 

THIS  PROGRAM  REGUI RES  A   SCRATCH  TAPE  FOR  ASSIGNING  CODE  VALUES  TO  LAYER 

7 

c 

DATA 

8 

DIMENSION  A1<35).A2(15).A3(10).A4(15>»STRAT(20.8)»A5(15) 

9 

c 

**THE  FOLLOWING  CARD  READS    IN  THE   NUMBER  OF  MATERIAL   TYPES    IN  EACH  RANGE* 

10 

RE AD (5.  2 00) I  NO » J NO , KNO » LNO » MNO 

1 1 

200 

FORMATI5I5) 

12 

c 

-***THE  FOLLOWING  CARD  PUTS  RANGE    1   MATERIAL   INTO  AN  ARRAY  FOR  LOADING*** 

1 3 

READI5.101)  (All  I  )  »I  =  l»INO) 

1  u. 

r 

***THE   FOLLOWING   CARD  PUTS  RANGE   2  MATERIAL   INTO  AN  ARRAY  FOR  LOADING*** 

1  5 
i  ? 

READ ( 5  1 101 )  (A2( I  )  .  I  =  l,JNO) 

1 6 

Q 

***THE  FOLLOWING  CARD  PUTS  RANGE   3  MATERIAL    INTO  AN  ARRAY  FOR  LOADING*** 

1  7 

READ  (5.101) (A3(I)»I=1  tKNO ) 

1  8 

c 

***THE  FOLLOWING  CARD  PUTS  RANGE   4  MATERIAL   INTO  AN  ARRAY  FOR  LOADING*** 

1  9 

READ ( 5 ♦ 101 ) ( A4 ( I ) #I=l»LNO) 

ZU 

c 

***THE  FOLLOWING  CARD  PUTS  RANGE   5  MATERIAL    INTO  AN  ARRAY  FOR  LOADING*** 

2  1 

READ ( 5  » 10 1 )  (A5( I  )  . 1  =  1 » MNO) 

22 

101 

FORMAT ( 10F8.0 ) 

2  3 

c 

***THE  FOLLOWING  CARD  READS   THE   TEST  HOLE  NO.i   NO.   OF  LAYERS  AND  MATERIAL 

c 

TYPE*** 

i$ 

10 

READ ( 8  » END =1001 ) NWELL .NOLAY  ,  STRAT 

c  o 

DO  1000  1=1.20 

?  7 

IF    ( NOL A Y . EG • 0 ) GO  TO  1002 

£  o 

r 

***THIS  SERIES  OF   TESTS  MATCHES  MATERIAL   TYPE   TO  ThCSE    IN  EACH  RANGE*** 

£  7 

r 

***WHEN   MATERIAL   TYPES  MATCH   THE   RANGE    IS   ASSIGNED  AND  STORED  ON  TAPE*** 

~x  n 

DO  20  J=l.INO 

•a  i 

I F ( STRAT ( 1 1 3 ) -Al t  J ) )  20.60.20 

20 

CONTINUE 

DO  30  J=l»JNO 

1  Zi 

It 

IF (STRAT ( I .3 )-A2 (J ) )  30.70.30 

3  3 

30 

CONTINUE 

J  fa 

DO  40  J=l.KNO 

1 

D  1 

IFISTRAT1 I «3)-A3( J) )  40.80.40 

3  ft 

40 

CONTINUE 

■a  o 

DO  50  J«=l.LNO 

/.  n 
H  u 

IF (STRATI  I » 3 ) — A4 ( J )  ) 5  0.90  ,50 

4 1 

50 

CONTINUE 

4  5 

DO  500  J=l»MNO 

"r  3 

I  F ( STRAT (  I  .3  1-A5  1J)  15  00  .5  01  .500 

44 

500 

CONTINUE 

45 

GO  TO  1000 

46 

c 

CODE    1.0.    2.0»    3.0.   4.0.   OR   5.0    IS  ASSIGNED  TO 

47 

c 

RANGE   1.0.    2.0,    3.0.   4.0.    OR   5.0  RESPECTIVELY 

48 

c 

RANGE  CRITERIA    IS  DETAILED    IN  TEXT 

49 

60 

CODE=1.0 

50 

GO  TO  100 

5  1 

70 

CODE=2.0 

GO  TO  100 

53 

80 

CODE=3.0 

a  i. 
54 

GO  TO  100 

90 

CODE=4.0 

5  O 

GO  TO  100 

C  "7 

501 

CODE=5.0 

58 

100 

DO  110  K=6,8 

59 

110 

STRAT (  I  ,K  )  =C0DE 

60 

1000 

CONTINUE 

61 

1002 

WRITE  i  9 )NWELL. NO LAY. STRAT 

62 

GO  TO  10 

63 

1001 

REWIND  8 

64 

ENDFILE  9 
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65  REWIND  9 

66  1004  READ(9»END*1003)NWELL»NOLAY.STRAT 

67  WRITE! 8 ) NWELL « NOLAY • STRAT 

68  GO  TO  1004 

69  1003  REWIND  8 

70  WRITE(6tl005) 

71  C     *»*THE  FOLLOWING  STATEMENT   INDICATES  RANGES  HAVE  BEEN  STORED  ON  TAPE*** 

72  1005  FORMAT (' ODATA  ON  TAPE  *  ) 

73  STOP 

74  END 

75  C  END  APPENDIX  D 


APPENDIX   E.-FLOW   CHARTS   AND  LIST  OF  VARIABLES  FOR 
PREPARING  GRAPHICAL  PRESENTATIONS 

Flow  Charts 


MAP  VIEW 


CROSS  SECTION 


READ  IN  TEST  HOLES  TO  BE 
PLOTTED 

t 

READ  SCALE  FOR  THE  PLOT 

XSCAL.YSCAL 

(SEE  VARIABLE  UST) 

1 

' 

READ  SPECIFIED  BOUNDARY 
OF  MAPPED  AREA  FROM 
DISK  FILE 

r 

SCALE  FIELD  DATA  POINTS 
TO  PLOTTER  COORDINATES 
(SYSTEM  DERNED  BY  USERS) 

t 

PLOT    MAP  OUTLINE 

_J 

j 

SUBROUTINE  WELPT 
READ  NW,  ALON,  ALAT 
FROM  CARDS  OR  DISK  RLE 
(SEE  VARIABLE  LIST) 


1 

r 

SCALE  PARAMETERS 
XLENG.  YLENG .  GTEST,  BTEST 
(SEE  VARIABLE  LIST) 

r 

RETRIEVE  LAYER  VARIABLES 
AND  STORE  ON  DISK  ,  NWELL, 
THICK  ,  KVALU 
(SEE  VARIABLE  LIST) 

\ 

r 

RETRIEVE  TEST-HOLE 
VARIABLES  8  STORE  ON  DISK 
NWELL,  GSURF,  TDEPH,  MAX  LA 
(SEE  VARIABLE  LIST) 

i 

CALL  SCALE 
SUBROUTINE 

PLOT 

CALL  GRID  AXIS  SUBROUTINE 
AND  ANNOTATE  THE  GRID 


PLOT  INDIVIDUAL  TEST-HOLE 
LOG  AND  WATER  SURFACE 


Figure  E-l. — Flow  charts  for  preparing  map  views  and  cross  sections. 
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List  of  Variables 


Variable 
name 


XLENG 
YLENG 
GTEST 

BTEST 


Format 


Description  of  variable 


Variable 
name 


Format 


Description  of  variable 


NWELL 
IVAR 

THICK 
KVALU 


NWELL 
GSURF 


TDEPH 
MAXLA 


Boundary  Card  for  Cross  Section 

F10.0  Distance  in  inches  between 
test-hole-log  plots. 

F10.0  Length  in  inches  of  the  cross- 
section  plot. 

F10.0    Reference  ground-surface 

datum  from  which  logs  are 
measured  on  the  plot. 

F10.0    Lower  reference  bedrock 

datum  beneath  which  test- 
hole  logs  do  not  extend. 

Lithology  Card  for  Cross  Section 

14  Test-hole  numbers  of  logs 
used  in  cross  section. 

14  Layer  number  for  lithology 
of  test  holes  in  cross  sec- 
tions. 

F9.1  Thickness  of  lithologic  layers 
in  feet. 

F9.1  Hydraulic-coefficient  range 
(1-5)  of  the  lithology  in 
each  layer. 

Test-Hole  Card  for  Cross  Section 


18 


F12.3 


F12.3 
12 


Test-hole  number  for  each 
log  in  cross  section. 

Ground-surface  elevation 
(feet     above     mean  sea 
level)  for  the  top  of  each 
log. 

Total  depth  of  each  test  hole. 

Number  of  lithology  layers  in 
each  test  hole  in  cross  sec- 
tion. 


Ground-Water-Elevation  Card  for  Cross  Section 


HWELL 
HCODE 


NYR 

NMO 

NDY 

NTIM 

ELEV 


15 
12 


13 


13 


13 


15 


F8.2 


Test-hole  number. 

Code  indicating  the  method 
of  measurement  of  ground- 
water level  (1  =  continuous 
recorder,  2  =  manual  tape 
down) . 

Year  when  water  surface  was 
recorded. 

Month  when  water  surface 
was  recorded. 

Day  when  water  surface  was 
recorded. 

Time  (military)  when  water 
surface  was  recorded. 

Ground-water-surface  eleva- 
tion (feet  above  mean  sea 
level)  in  the  test  hole.  One 
card  per  test  hole. 


Location  Card  for  Map  Views 


Note:  four  maps  are  plotted:  distributive  map,  sub- 
datum  map,  hydraulic-coefficient  map,  isopachous 
map. 

NW  F12.3     Test-hole   number  for  each 

hole  plotted  in  map  view. 

ALON  F12.3     Longitude  in  degrees,  min- 

utes, seconds  of  test-hole 
locations. 

ALAT  F12.3     Latitude  in  degrees,  minutes, 

seconds  of  test-hole  loca- 
tions. 

Lithology  Card  for  Subdatum  Map 

NW  13       Test-hole    number    of  each 

test  hole  that  penetrates 
the  specific  datum  plane. 

ISUB  15       Elevation  of  the  datum  plane 

to  be  plotted  (feet  above 
mean  sea  level). 

IDEPH  14       Depth  below  ground  surface 

to  the  top  of  the  lithologic 
layer  in  which  specified 
datum  plane  lies. 

ELTOP  F8.2     Elevation  of  the  top  of  the 

lithology  layer  that  con- 
tains the  specified  datum 
plane  (feet  above  mean  sea 
level) . 

RANGE  F12.4    Range   of  hydraulic  coeffi- 

cients (1-5)  at  the  speci- 
fied datum. 

THICK  F12.4     Thickness  in  feet  of  the  lith- 

ology layer  found  at  the 
specified  datum. 

DAPLT  F12.4     This   variable   name   is  as- 

signed to  the  range  of  K, 
S,  or  T  to  be  plotted. 

Test-Hole  Card  for  Map  View 

Note:  four  maps  are  plotted:  distributive  map,  sub- 
datum  map,  hydraulic-coefficient  map,  isopachous 
map. 

NW  15       Test-hole  number. 

THVAR  F5.0      Variable    THVAR    may  be 

any  of  the  21  test-hole 
variables  used  for  plotting. 

Net-Thickness  Card  for  Isopachous  Map 

WELL  NO.  A10      Descriptive  literal  field. 

NW  16       Test-hole  number. 

NET  THICK        A10      Descriptive  literal  field. 

THICK  F13.0     Net    thickness    in    feet  of 

lithologic  layers  having  the 
selected  hydraulic-coeffi- 
cient ranges. 
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Variable 
name 


Format 


Description  of  variable 


Net-thickness  Card  for  Isopachous  Map — Continued 

CODES'  A6       Descriptive  literal  field. 

CODE  (I)  5F4.0     Number  (1-5)  of  the  ranges 

of  hydraulic  coefficients, 
that  is,  K,  S,  or  T,  to  be 
plotted. 

Card  for  Hydraulic-Coefficient  Map 
NW  14       Test-hole  number. 


Variable, 
name 


Format 


Description  of  variable 


Card  for  Hydraulic-Coefficient  Map — Continued 


AKVAL 


F5.0 


ASVAL 
ATVAL 


F5.0 


F11.0 


Permeability  averaged  for  all 
layers  in  the  test  hole,  gpd/ 
ft-  of  aquifer.  The  size  of 
the  fields  containing  hy- 
draulic coefficients  may  be 
altered  to  accommodate  re- 
trieved data. 

Storage  coefficient,  dimen- 
sionless  value. 

Transmissibility,  gpd/ft  of 
aquifer. 
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APPENDIX  F.-INPUT  CARDS  FOR  GRAPHICAL  PRESENTATIONS 
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Ranges  1.,  2.,  3.,  1,,  and 
5.  ape  used  in  this  example 

TZ-89 
(5)3000 

X 

Z9-M 
(tj)3Q03 

X 

£9-09 
(9  3003 

x 

6S-9S 
(2)3003 

X 

55-25 
(T)  3(100 

X 

TS-9b 

X 

Sh-Obg 

o 

X 

Si-iZ 
MIHi 

93~ZIi — 

9T-TT 
MN 

CD 

OI-I  j 

X 

Nwmoo 

CROSS  SECTION 

DISTRIBUTIVE  MAP 

Q_ 
<C 
;  " 

e: 
«e 

ca 

HYDR.-COEF.  MAP 

ISOPACHOUS  MAP 

CQ  / 

<x.  / 

CC  /  O 

5/  UJ  — 1 

SZ-ST 
TVAiV 


W-0T 
1VASV 


6-5 


fj-T 


Nwmoo 


a; 
-> 

s 

60 
C 

o 

"E 
o 


9£-S2 
JATIV 

aaniun 
(oaoisi 


«-£I 
N01V 

3anii9Noi 
(waoisi 


3T-T 
HN 

310H  1S3H 
(T)aOlSI 


Nwmoo 


C3I —  — ■  : 

3^33; 


95-5h 

ndva 


ht7-££ 
>OIHJ. 


39NM1 


02-£T 
dOJ.13 


ZT-6 
Hd3QI 


8-tj 

ansi 


Nwmoo 


c i95o?si 


5-T 
310H  1S31 
(T)dOlSI 


Nwmoo 


T3 

u 


ft 
o 


ce:  cq 
<Ld« 


504  2  70  S  14  1346  1020.  10  "  

A  B  C  0  E  F  6  GROUND-WATER- 
  .    .   ELEVATION  CARD 


504.000       1037.300         43.000  5.000 
ABC 


504     I         8.00  2.00 
A     B  C  0 


TEST-HOLE  CARD 


BLANK  CARD 


LITHOLOGY  CARD 


'  I.  20.  1080.  940. 

A  B  C  D 

BOUNDARY  CARD 

Figure  F-3. — Input-card  setup  for  plotting  cross  sections. 


504  43. 
A  B 


TEST-HOLE  DATA  CARD 


504.000   974839.000  345751.000 
ABC 


LOCATION  CARD 


Figure  F-4. — Input-card  setup  for  plotting  test-hole  data  (see  figure  10). 


519  1080   70  1096.50 
ABC  D 


4.0000 
E 


23.0000 
F 


4.0 
G 


LITHOLOGY  DATA  CARD 


519.000   981205.000  350327.000 
A  B  C 


LOCATION  CARD 


Figure  F-5. — Input-card  setup  for  plotting  data  at  a  specified  subdatum 

(see  figure  12) . 


504  293.  0.003 
ABC 


HYDRAULIC- 
COEFFICIENT  CARD 


504.000    974839.000  345751.000 
ABC 


LOCATION  CARD 


Figure  F-6. — Input-card  setup  for  plotting  average  permeability  coefficients 

(see  figure  14) . 


WELL  NO.       504  NET  THICK 
A  B  C 


24.  CODES  T     T.    J^^.  £ 

D         E  F     G     H      I  J 

NET— THICKNESS  CARD 


504.000   974839.000  345751.000 
ABC 


LOCATION  CARD 


Figure  F-7. — Input  card  setup  for  isopachous  maps  (see  figure  16). 


TRADE  NAMES  ARE  USED  IN  THIS  PUBLICATION  SOLELY  FOR  THE  PURPOSE  OF  PROVIDING  SPE- 
CIFIC INFORMATION  MENTION  OF  A  TRADE  NAME  DOES  NOT  CONSTITUTE  A  GUARANTEE  OR 
WARRANTY  OF  THE  PRODUCT  BY  THE  U.S.  DEPARTMENT  OF  AGRICULTURE  OR  AN  ENDORSE- 
MENT BY  THE  DEPARTMENT  OVER  OTHER   PRODUCTS  NOT  MENTIONED 


*1976-G . P .0 . -1750-S/771-040/5 


.  S.  DEPARTMENT  OF  AGRICULTURE 

AGRICULTURAL  RESEARCH  SERVICE 
SOUTHERN  REGION 
P.  O.  BOX  5  3326 
NEW  ORLEANS,  LOUISIANA  70153 


OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE,  S300 


POSTAGE  AND  FEES  PAID 
U.  S.  DEPARTMENT  OF 
AGRICULTURE 
AGR  101 


